! 
\ 


Les solutions que nous venons d'obtenir permettent 


de résoudre les problèmes. électrostatiques qui: 


oncernent le système constitué par une grille et 
les deux plateaux symétriquement placés par rapport 
à celle-ci. Pour distinguer les uns des autres ces 
trois conducteurs, nous ferons correspondre l'indice o 
au plateau CD, l'indice 1 à la grille et l’indice 2 au 
plateau AB. Nous prendrons pour origine des poten- 
tiels celui de CD, le potentiel de la it étant U, et 
celui de AB étant U,. 


“2 Le problème que nous avons résolu à propos du 
. mouvement du gaz correspond au cas où U, est nul 
et où v, est remplacé par U,. Le second problème, 
celui de la distribution du champ électrique corres- 
pondait à U, égal à 2 U.. 


En superposant les solutions correspondant à ces 
deux problèmes dans une proportion convenable, 

on aura Ja solution du problème général pour notre 
système lorsque les potentiels U, et U, sont quel- 
_ conques. Si en effet U' représente la solution du 
: premier problème et U/ celle du deuxième, on aura 
évidemment 


NK 


(83) 


: LUE Ur; 


() Journal de Physique et le Radium, 1949, 40, 137 et 225. 
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SUR UN ANALYSEUR DE MOBILITÉ POUR LES IONS GAZEUX 
- III. APPLICATION DES SOLUTIONS DES PROBLÈMES TRAITÉS 
| A L'ÉTUDE DU CHAMP ÉLECTROSTATIQUE ENTRE UNE GRILLE ET DES PLATEAUX 


Par P. LANGEVIN. 


Sommaire. — Dans ce Membre: suite de ceux parus dans les deux numéros précédents (1), P. Lan- 
‘gevin calcule le champ électrique produit par une grille conductrice placée entre deux plans conduc- 
teurs parallèles. Après avoir étudié divers cas particulier il donne la solution pour une grille épaisse 
placée entre deux plans conducteurs situés à des distances différentes. 


d’où U égal à U,, et sur le plateau AB on a 
C0: TESTS 


d' où U égal à à Un. 
Les quantités les plus intéressantes à calculer 


parmi celles qui caractérisent électrostatiquement 


notre système sont d’une part la capacité C, de la 
grille par unité de sa surface et le coefficient d'in- 
fluence Co où C>9 des plateaux AB et CD l’un sur 
Pautre à travers la grille, Les coefficients d’influenceC;» 
ou GC de # grille sur les plateaux sont tous deux 


ÉgaUx à — — ; Par raison de symétrie. 


Grille mince, — Prenons d’abord le cas de la 
grille infiniment mince et de la distance d suffi- 
samment grande par rapport à la période a + b, cas 
auquel s’applique notre première solution du premier 
problème. La formule (30) donne 


= ui > : ) (2). 
d — D 


D'autre part, nous avons déjà remarqué que pour 
U,—2 U,;, la présence des fentes dans la grille 
infiniment mince ne modifie pas la répartition des 
potentiels, le plan de la grille prenant spontanément 
le potentiel moitié de celui du plateau CD. On a 


(2) L, logarithme népérien. 
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donc, en conservant l’origine des coordonnées ou .Y 
sur la grille comme dans la figure 5 (©. 


TRE m PAU > 
U = Us :) 


D'où, en portant dans (83) : 
eue += U; : (84) 


> Jo—Ly 2 


UV: (us 


Pour calculer les capacités ou coefficients d’in- 
fluence, il faut se reporter à la distribution des 
charges qui est déterminée par la fonction de flux © 
conjuguée de U. La fonction conjuguée de Y étant X, 
donnée par (24), et celle de n étant £, la fonction o, 
supposée s’annuler sur la ligne de force qui passe 
par l’origine, milieu d’un barreau, devient d’après (84) 


: ASUS SE UE 
Res 


1 


(85) 


Pour trouver les charges présentes sur le plateau CD 
par période a + b de la grille, nous devons donner 
à n la valeur à et à £ la valeur +, on a, dans ces 
conditions 


Æ=Ët=7, 


à étant supposé assez grand pour que Coth à puisse _ 


être considérée comme égale à l’unité. L’expres- 
sion (85) devient 


Le 
> }J d —LY 20 


D'après le théorème de Gauss, la charge Q corres- 


+; donc la charge 


LL 
7 


pondant au flux 9 est égale à 


portée par le plateau CD par période a + b et par 
unité de longueur des barreaux est 


7e 2 Üi— Us Un 
ECO LT) 17-806 


Par unité de surface de la grille, on aura la charge 


2 Ui — U b U; 


T#[rd— (a 6)L;y]  8rd. (5) 


SNA TE D 

Si l’on veut avoir la capacité de la grille, il faut 
supposer U, égal à zéro et doubler la charge Q, pour 
tenir compte de celle que porte en même temps 
le plateau AB. En divisant par le potentiel U, pour 
avoir la capacité, on obtient 


I 


ie 2[rd—(a+b)Ly] 


(81) 


ou, pour une surface S de la grille, et en remplaçant y 
par sa définition (22) 
r b : 
ed en 


Quand a est nul, c’est-à-dire que la grille est rem- 


S 


RES 
[sd (a + D) L sin - 


(:) Voir Journal de Physique Radium, 1949, 10, 186. 
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(A 


placée par un plan conducteur tete Je Rras 
est nul et l’on retrouve bien la capacité PE. d’un 


plateau isolé entre deux plateaux symétriquement | 
placés par rapport à lui à distance d. 4 
Lorsque, au contraire, b s’annule le logarithme | 
devient infini négatif et la. capacité s’annule comme 
elle doit le faire puisqu'il n'existe plus rien entre les 
deux plateaux. “4 
Pour trouver le coefficient d'influence du pla- = 
teau AB sur le plateau CD, par unité de surface de 
la grille, nous devons supposer nul le potentiel de 
celle-ci et (86) devient . 


ne Us I : 15 

Q=+ | sol 

d’où ; 
3 1 APE 

de Fe - à ce 
Ce coefficient s’annule bien quand la grille est pleine 
(a — o) puisqu'elle fait écran à l’influence mutuelle 
des plateaux, et devient bien, quand la grille dis- 
paraît (b — 0) égale à Ve 


I 
GES 8rd’ 


représentant l'influence, par unité de surface, : 
d’un plateau sur un plateau Pos situé à dis= a 
tance 2 d. a 


Grille épaisse. — os maintenant ce que 
donne notre seconde solution, qui est plus générale 
puisqu elle correspond à une distance d gere : 
et à une épaisseur quelconque des barreaux. 

il résulte de (40) 


= A5 + 2 +>4 res 2 cos ne, 
la fonction conjuguée : 
tu SA CRI Rare 
PE sx LEE sm27n6, 
et de (67) - 
U"= Bo à +S 2, Re ee - 
ù : sh2nû 


la fonction con) ee 


o" LR >: a ; 6 
= ee sin2 RE. 
: ee 


La charge du plateau CD par période a + b et par 
unité de jORBUeEE des barreaux s’obtient en faisant 


dans les w, £ égal à x, d’où, par application de la 
formule déduite de (83) 


la charge 


Fe  - 
List, f 


: re 


et, F 


RE o U; ” A5 PU B5 
dE 4rd 2 U; Ard 


; _ Le cas où la large des fentes est égale à la lar- 
Eur des barreaux se présente comme particuliè- 
rement simple, puisqu'il suffit, pour avoir A, de 
résoudre le système d’équation . constitué par la 
première équation du groupe (55) et les équations 

du groupe (56), B, s’obtenant de la même manière 
à partir de la première équation du groupe (72) 
4 et l’ensemble du groupe (74). 


capacité C; et pour le coefficient d’ HORERUE à tra- 
vers la grille C,, 


A5 et C a A9 B, 
2rdU SUR: 


CE 


Nous trouvons là encore une occasion de mettre à 
: j épreuve notre méthode d'approximations succes- 
Be en confrontant les résultats (90) avec ceux de 
+ première méthode qui deviennent, dans le cas 

L rce eee ici où .a est égal à b 

Ce 

4 —= HER RD TENN ETS 
 2(rxd+0%al V2) 

F (90”) 


2 ch I D es 
_ SR F0 rd+oal Vo sd. 


= Pour ce qui concerne C;, l’approximation que nous 
avons déjà utilisée pour os les résultats (62) 


ES AE Ô 

LA Ô Æo 413 

‘où, d’après (90) 
Ci = 


te SU ,834a) 


gx 
7 au lieu de, d’ après (90) 
n- 
C1 = 


I 
et 6924) 


nul et le groupe (74), dont les équations deviennent 
homogènes en B;,,; avec un déterminant non nul, 
- donnent des valeurs nulles pour ces B,,,, de sorte 
À que la première équation (72) nous conduit à 


Bo 
U 


c I I 1 
D — - ) 
#5 8rd+0,83%a HA A 


Y I ! c 1 
Gun ere | 


L'approximation suivante, où l’on conserve trois 
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nu 


La formule générale (86) nous donne, pour la 


(90) 


259 
des Ap+as Soit A, A, et À;, donne 


LI 


C1 = 
1 ; + 0,7814&) 


et 
I ll Il 
C> = ee ee M 
Ë ne | 
La suivante encore, avec quatre inconnues d'’in- 


dice 1 impair, À, 43, À; et À dans les équations (56). 
et la première des équations (54) donne 


ll 


OS NS 
2(rd+0,7604a) 


et 


PL 1 Ù 
CE k 
e ere | 


On voit ainsi les approximations successives 


converger vers la valeur exacte. 


À 


Problème électrostatique général de la grille 


entre plateaux parallèles. -— Nous nous sommes 
limités jusqu'ici au cas du système symétrique où 


les plateaux sont placés de part et d'autre de la 
grille à égale distance de celle-ci et la solution (83) 
pour des potentiels quelconques des trois conduc- 
teurs a été obtenue par superposition de deux solu- 
tions particulières U' et U” qui correspondent aux 
deux cas que nous pouvons appeler respectivement 
le cas symétrique où les plateaux sont tous deux au 
même potentiel pris pour origine, et le cas anti- 


symétrique où les différences de potentiel entre la 


grille et les deux plateaux sont égales et de signes 
contraires. Il nous sera utile d’expliciter ces deux 


. solutions particulières U! et U” en fonction def et n 
P ; 


ainsi que les fonctions conjuguées 9’ et 9’ qui leur 
correspondent. Si nous conservons les mêmes axes 
de coordonnées (fig. 6) et les mêmes notations que 
précédemment, à ceci près que, dans la solution 
symétrique la vitesse v, de la grille sera remplacée 
par son potentiel U;, la fonction U' est donnée 
entre CD et GH, c’est-à-dire pour nñ variant entre o 
et 9, par l'expression analogue à (40) 


DE sh2r sh2n" 
UE AA 5 me ñ cos 21e. 


Q 
sh2nÔ 2e) 
La solution entre GH et AB, valable quand x varie 
entre à et 2 9 doit être symétrique de la précédente 
par rapport à GH; on a donc, dans l'intervalle 
considéré 
= AS ee +Y 4, RÉ cos2n€t. (9) 


sh 2 7 à 


Cette solution s’annule bien sur AB, et les deux 
fonctions (91) et (92) se raccordent bien sur les 
segments de GH situés à l’intérieur des fentes sans 
présenter aucune HUE dans les valeurs de U 


et dans celles de = qui y est nul. 
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La fonction 9 étant relié à U par les conditions (75), 
on voit qu'il correspond à (91), pour la région com- 
prise entre CD et GH, c’est-à-dire pour o <n < Ô 


= 1e +Y A4, 


et à (92) pour la région 9 < n < 20 


a h°2 DO) 

Ë a ch2n(20 fra 
= NA EE sinon E. 
< (9) 


sh 2710 


ch2nn 


(93) 


sir 2 n E 
sh2 n d 


(9%) 


De même, la solution antisymétrique donnera les 
formules correspondantes 


l Eee sh2r eee 
Bert 2 PB 50 OS27NE, (95) 
Ô sh 2 7 Ô 
Bo(28 —") 
[#1 = or = 
Q 
sh2n(20 — n) F 
Es © Ph AS PA D COS ONDES (96) 
sh 2 nr à es 
ë ch2n"n . E 
= BE +Ù 8 ” ll énone, C9) 
! ù sh2 nr 0 
£ h2n(20 — nn) 
12 GE Y cn? \ l RE aR\ 
D Bo = + B,, ————— sin2neE. 98) 
; PSE dl > sh 2 7 û É ( 


Ainsi que nous l’avons vu précédemment, la 
fonction de flux © permet de calculer les charges 
présentes sur les plateaux; pour le plateau CD, 
il faut faire n égal à o dans les formules (93) et (97) 
et prendre la valeur de © pour £ —rT; en désignant 
par Q, et Q% les charges présentes par unité de 
surface de CD dans la solution symétrique et la 
solution antisymétrique respectivement, on obtient, 
en grandeur et en signe 


0! dE o PES A 
USE ete eh) -{rd 
; (99) 
De Pr 
si0 ir(a+b) Ir 


On aura de même sur le plateau AB, en remarquant 

que la normale extérieure y est dirigée dans le sens 

négatif des n 
QUE 


#1) 4 
Nr À 


Ax(a< bd) 0. And. ‘ 
: ) (100) 
C PRES Te OF RT B f 
SAT Er(a+b) 4rd. 


La charge de la grille par unité de surface est égale 
et opposée à-la somme des charges présentes sur les 
plateaux, donc 


Pa 
Te 


Lorsque la grille et le plateau AB ont des potentiels 
quelconques U, et U, par rapport à CD nous avons 


# 
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vu que la distribution du potentiel est donnée pe ir (83) 


% R à 
LU DEEE U: 7 H 2 W 
= fs 7) v cu. 3 
‘ On aura donc ? 
| . Us RE Us 1 
— = GERS U; - = U; 
et pour les charges, en posant = 
A) PB | 
= — bo= — 101 
do T: ? 0 U; } ( L ) 
on obtient 4 
SA do Do 5. do 
EP PE 4rd Ds 
à 07 bo ee 0 ; f: / n9 , 
Ga at 
(en y 
= 274 Rd: 


Il en résulte en particulier pour la capacité de 


la grille et pour le coefficient d'influence mutuel des 


plateaux à travers la grille les valeurs, par unité 


de surface. 


Les coefficients a, et b,, purement numériques, 
s’obtiennent par la résolution des systèmes d’équa- 


tions linéaires établis pour la solution de chacun des 


deux problèmes symétrique (aç)etantisymétrique(b;). - 


Quand d est suffisamment grand, il résulte de la 


forme particulière que prennent alors nos groupes. 


d'équations linéaires, que a, et b4 se présentent sous. 
les formes : 


ù b Ô 
— Æ ———— 
Ô+a 570 BY 


y —= 


où les quantités « et 6 dépendent uniquement de 


la grille, par l'intermédiaire de & et de o. Quand «. 
est nul (grille infiniment mince), nous avons vu 


que les valeurs limites de « et 5 sont 


À Fr 
&=Téin 0; ë 6 0: 
2 


Quand 9 est nul, la grille disparaît, « devient 
infini et 5 nul, de sorte que a, est nul et b, égal à 
égal à 1, les fentes disparaissent, la 


l’unité. Pour p 
grille devient se * plaque continue, « et 5 deviens 
égaux à —6. 

En substituant la valeur (104) de a, dans l’ex- 


pression (103) de la capacité C;, on obtient pour 


celle-ci 
I 


ER 
: a(a+b)] 


5[rd— (105) 
avec ce résultat que l'influence des fentes sur la 
capacité de la grille, équivaut à une augmentation 
de la distance des plateaux de chaque côté de la 


403) 


06) 


a (ab), 


74 La simplicité de ce résultat est liée au fait que, 
4 lorsque la distance d devient suffisamment grande 
par rapport à a + b, le champ devient uniforme 


: tions, pour une même distribution du champ au 
voisinage de la grille, et par conséquent pour une 
même distribution des charges à la surface de celle-ci, 
… un accroissement de la distance et produit simple- 
: ment un accroissement proportionnel dans les 
_ différences de potentiel U, et U,. Ces différences 
. de potentiel sont donc des fonctions linéaires de la 
distance quand celle-ci devient suffisante, et c’est 
bien ce que traduit, en particulier, la forme (105) 
_ dela capacité de la grille. 


Cette uniformité du champ au voisinage des pla- 

teaux permet également de résoudre le problème 

de la distribution du potentiel pour des valeurs 

. données de celui-ci sur la grille et les plateaux 
porte ceux-ci sont à des distances différentes de 
_ la grille, ces distances d et d' restant suffisamment 
| grandes par rapport à à la période a + b. 


_ Soit A’ B” (fig. 6) la position du plateau supérieur, 
_ à distance d’ du plan médian de la grille, le plateau 
- ténor restant à distance d. Soit U’, le potentiel 
_ sur le plan A’B', potentiel qui sera constant dans 
. toute l’étendue de ce plan, en raison de l’uniformité 
du champ entre AB et A'B'. Ce champ a pour valeur, 
en vertu de (83) 


2 ( U; Le User 
D SU) 07. EU O7 


- (406) 


-_ Les valeurs à prendre pour U' et U'!' sont celles 
__ qui conviennent pour la région supérieure du 
* domaine, et qui sont données par les formules (92) 
et (96) respectivement. On en déduit 


aU" A5 DUC: k; 
dy a dy d? NES 
_ pour les valeurs de n voisines de 20, où shond est 
grand. 


De même on déduit de: (106), en y portant les 
_ valeurs (107) 


oU 


5 (D 
i + == Ca bo) — RTS 


d 


On a d’ailleurs, puisque ce champ est constant 

entre AB et A’B' 
/ OC 
UE ed a Ve 
CU 


L4 — 


el 2) + Tata) 
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au voisinage des plateaux et que, dans ces condi- 


d’où 
: ad U,, RS (= 4") 0 U 
2d— (a+ bo)(d— d') 


U; — 


Il suffit de porter cette expression de U, dans la 
formule (83) pour avoir la distribution du potentiel 
entre la grille et les deux plateaux A’B' et CD 
lorsque le potentiel de A'B' est U’, et celui de la 
grille U, celui de CD étant toujours pris pour 
origine 

dUs—(d— djaoLi UV'—U 


(108) 


us rt 
U C ae 4 — (d — d") (&i 5. bo) U, 
La fonction de flux o sera par conséquent 
TUE d'te Un oo à 
FE te de ant) 0 


En tenant compte de (99) et (100), on obtient 
pour les charges par unité de surface des plateaux 
et de la grille 


f(ao— bo)dU ) 
De lof ad = bit d = UN 
Van dr de (are bdd) 
(as + 60) US = 2&0 Ui 
Gr[2d = (a b}(d =) 
[2404 — avbo(d — dd) U— &dU, 
2rA[24 — (ao + by)(d — d')] 


) 


(10) 


Oo 


= 


De là résultent les valeurs de la capacité C, de la 
grille et du coefficient d'influence mutuelle des pla- 


teaux à travers la grille Cy 


244 = 4 by(d — d') 
2rd[2d— (a+ b)(d—d)] 
y — bo 


inl2d (a+ bdd) 


Ci 


(114) 


C5 


Quand d est égal à d' on retrouve bien les expres- 
sions (103) que nous avions obtenues dans le cas 
symétrique. 


Enfin, dans le cas de la plaque pleine, lorsque . 
fentes de la grille ont une largeur « nulle, les cons- 
tantes « et 5 des formules (104) sont égales à — e, 
on a donc 

Ô d 


(44 UN TAPIE 
9 UNE RIÈESTE 


et les formules (111) deviennent 
I I 


Ca rdc) — Fed c) 


Co + 0, 


comme il est raisonnable, puisque la plaque est à des 
distances d — c et d' —c des deux plateaux, et que 
cette plaque pleine fait complétement écran à l’in- 
fluence mutuelle des plateaux situés de part et 
d’autre d’elle. 
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Traitement direct du cas dissymétrique par 
la méthode des séries trigonométriques. — Je 
vais, 
méthode des séries trigonométriques le cas général 
dissymétrique de la grille placée entre deux plateaux 
situés à des distances d et d' de part et d'autre du 
plan médian de celle-ci. Soient de nouveau o, U;, U;’ 
les potentiels imposés respectivement au plateau CD, 
à la grille et au plateau A’B’. La distribution du 
potentiel entre CD et le plan médian GHsera donnée 


par ; 
WA 0 . + A 
(e] 


La distribution entre GH et A’B' sera 


S/ 
d 0e " 
CPE A" 0 ii 


Ô 

., Shon(d+5—"n) 

— PRET Ton RE 
sh2nû 


u s’annule sur CD en même temps que n et u 
devient égal à U, sur A'B' pour n égal à 9 —0!. 
Les conditions à remplir par les 4, et les A’, doivent 
traduire les faits suivants : 


sh2n1n 


n = COS? 7€. 
sh n 0 < 


| (119) 


cos2nE. 


! 


10 u doit être égal à U, sur RS en même temps 
ue u! est égal à U, sur VT:; 
1 ; 


20 u doit être égal à U, sur SN et u' FES à U, 
sur NT; ; 


30 Les deux fonctions u et u! doivent se rac- 
corder sur NM' de façon que u et u' y prennent la 


même valeur en chaque point ainsi que leurs dérivées 


premières par rapport à n et à £. 

Ainsi les deux fonctions harmoniques se raccor- 
deront puisqu'il n’existera sur les bandes médianes 
des fentes, telles que NM’, aucune couche double 
ni simple correspondant aux surfaces de disconti- 
nuité du potentiel ou du champ. 


Nous aurons ainsi trois groupes d'équations 
linéaires auxquelles doivent satisfaire les A, et les 4’, 


I. Le premier groupe sera analogue à (43) et 
s’écrira 


à Ô —€ I An Shon(d — e) : 
CNE AN SE > 2 ‘ —— Sin À TE, 
Ô To n sh n ù d 
nA} shont(on ru 
Ÿ sin 270 — 0: 
Zn 127 — mn? sh2nû Ë e 
Rs 1 
s d'— + I A c 
Bet A; — ù —" (113) 
Ô TP 1) PIN 
sh 27 (0%): 
a. OI 70 70 
sh 7 0 se 
X n A}, sh3n(d +) - 
2. Sin 2 70 = 0. 
Le N°2 0? — n° sh2nû a 


IT. La seconde condition conduit à des équations 


pour terminer, traiter directement, par la. 


analogues à (51) 


ch2nû-ch2n(ù—e) ; n 
COS 71 T9, d 


sh2n0ô ©: 
RS €) 
sh 2 7 à 
=-X COS TE.— 0, 


sh2n(ù — 2] 


sh 2 7 à 


= 
Le A n Pis 
di N°2? + m7? 
NA ; C 

XX COS RTE < 0: 


ch2nd — chon(à — ec) 


Shane COS 7 TP, 


C2 


4 TA 
= nA,, 
>. n?e?+ mr? 


l 
ge 1 | 
É sh27n(à —e) 
ER mer 
Teil sit 2720; 


ch2n0 — ch2n(ù — :) 
sh2nû 


! 
D A 
n2e? + mn? 7? 


l 


X COS2TO = O0. 


III. La continuité entre u et u' le long de NM 
donne 

NES / : I ; < | sn 

CA ms : d $ 


A 


> Cned 
. Le Fe = =. 


os ds. 


(LISE 
s ANTITSS PER : *e - 


FRS . du 
La continuité entre — 


ou 
D. et FA le long de us 


donne : | 
A5 ; f 1 : | : | ; 
Bras) > 1) 
x (A, coth2 70 — A4;,coth on 5) = 
KSsnnT(r — = ‘& 
FAT (116). 
NE | 


di °2(1— 9) — me? 
x (Ancoth2nù — coth2nù') 


& X SiNnT(I—p)= Oo. 
du 


=. le long de NM' 


La continuité entre 


donne 


D 


R?(1— p }} — m? 


(117) 
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En liaison avec la distribution (112) du poten- 
tel, celle de la fonction conjuguée o prend dans les 


deux régions situées de part et d’autre du plan 
- médian se les valeurs 


. chon 
DAS ñ Ta sin onE, 


< 


mn = A 


CAE NY 


< e ch2n(o 0 "1: 
3 = 45 + DA, chon(5 +5 —"n) 


ons Sin 2 n €. 


_ teaux. cb et AB" 


0 A 
Qo= — Ar d Q:— rat 


_ La résolution de notre système d'équations 
linéaires donnera A, et A en fonction linéaire 
de U;et U, et l’on déduira immédiatement de (119) 
- les coefficients d'influence et les capacités. 


extrêmement simplifiée dans le cas où la largeur a 
des fentes est ei celle b des barreaux, c’est-à- 


dire où o est égal à: = 


JLest en effet Dole dans ce cas, en procédant 
comme nous l’avons fait antérieurement, de former 
un groupe d'équations linéaires ne renfermant plus, 
en dehors de A, êt A’, que les inconnues d’indice 
3 impair A1 et ae. Ce groupe est le suivant 


Fe +. ST PS 
4 (BrEE = A Ce ) L Asp, | 
EN 7 2P 5-1 
Ô > —1)?,S ; 
US At EE (( ) As 
cu 2 + TI 
e 0 
2 (—1} 
=. + A = À + À, 2e 
UE À 0 Mers p+H1 P= 1) 


_X[coth2(2p +1)0 Asp 
+ coth2(2p +1)0 A°p+#1], 
(ir (2p +0 : 
En de Le 


(120) 


> [ Se Smp)A2p+1 = 710 que Sp) A2p+1] 29; 


(—1}(2p +1} 
DOTE 


$ : om 
= X 
: È fée HT 


+ coth2(2p +03) A2p+i 


coth{môS np . 


2m SA 
— coth{ma Sy | 
2pP Ca m0 Le 
:] 
| 
| 


Ap | 


+ coth2(2p + 1?) 


à ee. SUR UN ANALYSEUR DE MOBILITÉ POUR LES IONS GAZEUX 


5e (118) 


__ Comme précédemment, on en déduit pour les. 
charges Q5 et @ par unité de surface des pla- : 


(19) 


La résolution du système d'équations se trouve 


les p RONA GT 
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Avec les notations 


Ke sh2(2p +1)( —e) sh 4 m û 
PO ISh 20 DE To sh 4mn(0 — +)? 
NE sh2(2p +1)(à —e)  sh4{mo 
A sh2(2p+1)0  sh4m(o — €} 
sh2(2p +1)(d—e) -ù 
Sop a = mo) 
sh2(2p +1)à Ô — € 
, sh2(2p +1)(5 —e) à’ 
Sop = RE 


sh2(2p +1) à —e 


L’élimination de A, et A’, entre les quatre pre- 
mières des équations (120) fournit les deux équa- 
tions 


(Sop # 1) A1p+1 


HS (Sp + 1) Aop+t], \ 


+ coth2(2p + 10) Apr 


Se 
d (2p +1) 
+coth2(2p + 1) 


*< Ars] . 


Ces dernières équations, jointes à celles des deux 
dernières lignes de (120) permettent de déter- 
miner les A4,,,, et les A°,.,, qui, portés dans les deux 
premières équations (120) donneront A, et A 


Une approximation déjà suffisante consistera à 
donner à p les valeurs o, 1, 2, et à prendre deux 
équations dans chacune des deux dernières lignes 


de (120) en donnant à m les valeurs r et 2. 


Pour aller plus loin et en choisissant pour p 
puce "qui. donne 
2(p + 1) inconnues 4:,:1 et A°,,,, On donnera à m 
les valeurs r, 2, …, p dans les deux dernières lignes . 
de (120), ce qui fournira avec les deux équations (121) 
les 2p +2 équations nécessaires pour déterminer 
ces 2 p + 2 inconnues. 


Travail fait à Troyes én 1942-1943. 
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SUR LES ÉQUATIONS D'ONDES DES CORPUSCULES DE SPIN 


PAPE 


SÉRIE VIII, TOME X, OCTOBRE 1949. 


ÂT 


ÉT LEURS SOLUTIONS DANS LES CHAMPS NUCLÉAIRES GÉNÉRAUX 


Par M. GéÉRARD PETIAU. 
Institut Henri Poincaré. : : RE 


ss 5 + 0 
Sommaire. —- Étude de l'équation d'ondes relativiste de Dirac représentant le corpuscule de spin - — 
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dans le cas où ce corpuseule est soumis aux champs de forces que l’on considère en théorie du noyau, 

plus généraux que le champ électromagnétique, et s’exerçant sur le corpuscule par l'intermédiaire 

de ses diverses grandeurs caractéristiques représentées par les seize opérateurs matrices de la théorie du 

corpuseule de Dirac. Après avoir discuté les approximations non relativistes correspondantes et notam- Fes 
ment les caractères des couplages spin-crbites, nous examinons en détail comment l'extension des TA 
méthodes de Dirac et de Darwin permet la résolution de l'équation d'ondes par séparation des variables 

dans le cas des interactions centrales. L'existence dans trois cas remarquables de l'intégrale première 

de Dirac permet par l'extension d’une méthode de Temple de calculer les niveaux d'énergie dans trois 

cas particuliers et d'écrire explicitement les solutions au moyen de combinaisons de fonctions radiales 

et de fonctions angulaires dans le cas d'interactions d’intensités constantes. 


I. — La représentation de l’électron par des 
fonctions d’ondes solutions de l’équation de Dirac 
à permis d’établir une théorie cohérente et en accord 
avec les résultats expérimentaux pour l'atome 
d'hydrogène, tout au moins lorsque l’on néglige 
les effets d’origine nucléaire, ce problème étant 
ramené à la recherche des solutions correspondant 
aux états d'énergie stationnaire d’une charge élec- 
trique élémentaire dans un champ électrostatique 
coulombien. 


Plus généralement, on admet aujourd’hui que 

l'équation d'ondes de Dirac permet de représenter 
/ 

tout corpuscule de spin et notamment le pro- 
1T 

ton, le neutron, et même certains types de mésons 


auxquels on attribue le spin _ De telles particules 


lorsqu'elles sont chargées subissent les actions élec- 
tromagnétiques de type classique s’exerçant par 
l'intermédiaire des potentiels électromagnétiques 


èt en outre des actions directes des champs électro- 


magnétiques par l'intermédiaire de leurs moments 
électromagnétiques propres. De plus, ces particules 
sont susceptibles de subir des actions d’un type 
plus général que nous appellerons nucléaire ou 
mésique. L'examen du comportement des parti- 
cules dans ces champs généraux et en particulier, 
la recherche des états stationnaires, nous paraît 
nécessaire pour l'édification d’une théorie des 
systèmes nucléaires et la discussion des conditions 
de la diffusion des particules élémentaires par les 
noyaux atomiques. 


; : 0 
Nous représenterons le corpuscule de spin ;—» 
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en l’absence d'interaction, par les fonctions d’ondes, 


solutions sous-les conditions ordinaires de la Méca- 


nique ondulatoire, de l’équation 
[po (px) + moca;]d =, (1) 
avec 


TOO) SC OPSE 
Dies Re PERS Fe 


ot 


d, (@ys os a), &4 FOTnent un eu de quatre 


matrices du quatrième rang telles que 
CE J =1,2, 3; 4) 


XX j + a ja; = 2 0j. 


que nous représenterons suivant la méthode de. 


Dirac au moyen de matrices du second rang o et p, 
soient 93, Te, Ta Pis Pas 0 telles que 


Spy = Sr PpPq—=—ÎÜr  (p;4; 


en posant 


Cahier (S hier  (mhjæi= Gi(ps);e 
(UT RAR) 


Nous adopterons pour les matrices g et p la repré- 
sentation 


OA 


61 ou P1 = 


Les matrices «, «, permettent de former un système 


À 


HR NP RNR EE S 


Me 16 matrices linéairement indépendantes repré- 
sentant. les grandeurs ER du corpus- 


_cule de spin e 


hs Xo—= 1], œ ; 


< . 
ü F CPR 
Go l4%a; = pi, 


PE = 703; 


(ep = ay); 


; AR TALE— TE 3P», Hs = 01 Lo ARR y = Le 


Dans le domaine où s exercent les interactions que 
- nous allons considérer, l’équation (1) sera rem- 
pce par. 


[po+(p.x) mes + H]$=:0, (4) 


CH, étant un opérateur d° interaction que nous écrirons 


Eine le cas général 


Hi= Las EU ( U.x)] 


L, est un invariant par rapport aux transformations 
de Lorentz: £ 

Us U un quadrivecteur: - 

 P, M un tenseur antisymétrique du second ordre; 

ES. S un pseudo-vecteur; 

a LE un pseudo-invariant. 


Nous nous proposons ici de discuter la possibilité 
d'extension dans le cas de l'opérateur d'interaction H, 
- de la théorie (1) que l’on développe dans le cas où 
2 cet opérateur se réduit au terme 


Ce problème a déjà été abordé dans des cas par- 
_ticuliers: ‘par de nombreux auteurs. C’est ainsi que 


l'introduction du terme (U.o), avec U = A, A poteri- 


tiel vecteur correspondant au champ dos 
dipolaire du noyau, a été considéré par G. Breit, 
dans l’étude de la structure hyperfine des niveaux 
de l’atome d'hydrogène (2). L'introduction du terme 
en (P, M), dont la possibilité a été signalée par 
* Pauli (®), a été utilisé par Bethe (#) pour étudier 
les possibilités d’ionisation d’une particule neutre 
(en l’occurence le neutrino) possédant un moment 
électromagnétique propre. Ce terme a été également 
considéré par G. Breit (5) dans la théorie de l’électron 
pour discuter l’effet du moment magnétique intrin- 
sèque introduit par J. Schwinger (6). Nous avons 
- discuté en détail la théorie de l’électron complétée 


- (1) P. A: M. DrrAc, Proceed. Roy. Soc., 1928, A 117, p.610; 
1928, À 118, p. 351. L. DE BroG1tE, L’électron magnétique. 
- (2)-G. BrerT, Some effects of the intrinsic magnetic moment 
- of the electron (Phys. Rev., 1948, 74, n° 6, 656-663). 
@) W. PaAuLzT, Handbuckh der Physik, 24, 1, p. 233. 
(*) H. BETHE, lonisation power of a neutrinos with moment 
magnetic (Proc. Cambridge Phil. Soc., 1935, 31, p. 108-115). 
(5) G. BrerT, Does the electron have än intrinsic magnetic 
moment (Phys. Rev., 19/47, 72, p. 984; 1948, 78, p. 1410-1471), 
(5) J. SCHWINGER, Phys. Rev. 1948, 73, p. 415. 


+IP:x) <CMp)l + [Soc +(S. 5) Bas (3) 


At : Er, 
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Dar une interaction électromagnétique de ce type 
dans un travail précédent (?). 

L'influence du terme en /,«,, dont la possibilité 
de l'introduction avait été signalée par Furry (8), 
a été discutée en détail par Pais (?), puis par Caldi- 
rola (1°) qui a déterminé at les niveaux 


d'énergie correspondant au cas 1, — = À étant une 


constante. 

Signalons également que, dans un travail précé- 
dent (1), nous avons discuté le problème du passage 
d'une onde corpusculaire solution de l'équation 
de Dirac à travers une barrière d'interaction générale 
correspondant à H,, avec toutes les grandeurs d’in- 
teractions Z,, U,, U, P, M, S,, S, I, constantes. 
Dans une Note aux Comptes rendus de l’Académie 


des Sciences (2) nous avons également discuté briè- 


vement l’équation d’ondes correspondant à l’ap- 
proximation non relativiste de l’équation (4). Nous 


préciserons ici les résultats de cette Note et nous. 


discuterons en détail le cas des interactions centrales. 


II. — Nous allons examiner en détail les équa- 
tions d’ondes correspondant à l’approximation non 
relativiste de l’équation (4). à 

Si nous supposons les interactions indépendantes 
du temps, posant 

lypr 
V2, 2 0)= (2 7: 2)e ; 


l’équation (4) s’écrira 


2 + (px) + Moca, 
ee Ar Ge, 
+ [(B:7) +(M.)] 


Lise bise loc 35) 07 06) 


Suivant la méthode de Darwin (#) nous décom- 
poserons ce no en deux systèmes de deux 
équations of! ou o,, 9? ou , en mettant en évidence 
par l'indice 1 ou 2 l'indice correspondant aux 
matrices 0. Se 


i e » 


(7) G. PerraAu, Sur quelques propriétés des corpuscules 
de spin = 2 dans les champs électromagnétiques (J. Phys. 


1948, 9, p. 228-24). 

(#) W. H. Furry, Phys. Rer., 1986, 50, p. 784. 

() Pais, Phys. Rev., 1945, 68, p. 227. 

(2) P. CALDprrRoLA, Sul moto di una particella di Dirac in 
un campo coulombiano misto (Nuovo Cimento, n° 1, février 
1948,.5, p. 1-7 

(2) G. PerrAU, Sur le passage des ondes corpusculaires de 


spin (2) (2) à travers les barrières de champs nucléaires 


(Revue Scientifique, 1947, 85, p. 1094-1r06). 
(2) G. PETIAU, Sur l'équation d’ondes non relativiste des 


/ Free ae 
corpuscules de spin = dans un champ nucléaire général 
LC: 
(C. R. Ac. Sc., Paris, 26 juil. 1948, 227, p. 263-264). 
(3) G. DARWIN, Proceed. Roy. Soc., 1928, A 118; p. 654. 
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Nous écrivons donc, au lieu de (6), 


[mme Dee h ets M). 1e 
+ [(œ+u —1P). e)+ Si |s Ps Die 1) 
14 

[5 — moe + U— + [(S — M). s1+: 
[ce HUE :Pj:s)= nr So ]gi=o. (8) 


Nous passerons à l’approximation non relativiste 
en mettant en iHeece in énergie non relativiste E 
en écrivant 

WE Ernie 


1 


BE . 
et en supposant : négligeable devant m,c. 


Si nous supposons en outre U,, 1, S et M nédgli- 
geables devant mc (hypothèse de faible couplage), 
l'équation (7) nous donne immédiatement 


Ll 


DIR — 


[Gæ-u-ipie) ess Do. (9) 


20 C 

et en reportant cette expression dans (8) nous 
obtenons pour 4, | équation 

F 


[l 


2 M0 C Fi U} + Si+13-+pP? 
= [(I2PAp+U) 


+ 2$0(p+—U) 21P).e|| 


_|aivp + [(RotU — :RotP 


20 C 


+1 grad Si— grad/:).6] | | @» —=0. (10) 


) 
Si nous introduisons l'impulsion généralisée 
H=p+U—(P/;)+So5, (11) 


l'équation (10) s'écrit encore 
E 
LE = + U— Li (cs = M).c) 


CE 


_— [divP + (RotU— gradr.).) || +2 =0. (12) 

Nous obtenons done, dans le cas de faible T. 
et à l’approximation non relativiste, une équation 
de Schroedinger-Pauli généralisée contenant un 
terme d’énergie potentielle dépendant du Spin, 


le terme d’énergie cinétique 


I ; 
H? dépendant 
270 € 


également du spin par l'intermédiaire d’un vecteur 
impulsion généralisé II dépendant explicitement 
du spin. 


Afin de préciser et d'approfondir ce résultat, 
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nous allons examiner Rues cas dre É 


3 
notamment les cas ou les fonctions d’ interaction | 
ne dépendent que du rayon vecteur (interactions 
centrales). ; : TEST 
19 Nous examinerons d’abord jé cas où 
Hi Tia + Uo+(U.x). 


AR 


L'équation (6) se réduit à ce Se 


3) 


[= + (p: QE Mo Ca + id + Vo + (U:. os 0; 


Si en particulier D:= grad #, “ÿ étant un inva- ‘ 
riant, ce qui sera notamment le cas si U est central” 


[u - U (r) — a pe 


U(æ, Dau de 


ee. 


l'équation (13) se raménera à 


[2 + Cp: &«) + MocCu Le +0 ]#= Un aa 


C 


Par suite, l'interaction vectorielle de ce type 4 
notamment l'interaction vectorielle centrale, peut 
être écartée par une modification du facteur de ? 
phase. 5 
A partir de l'équation (13) nous obtenons les 
deux équations en ,, % Fr 


|» Mo C + - lo n|e oi + [(p + U).318. 0 cie 
* [S+n-rer(p-urses Gé). 

Si E = mé, nous en tirons | | 
A 


2m0C + Uo + 


(nous ne faisons pasicil’hypothèse de faible couplage). 

Reportant cette valeur de v,, l’équation (6 À 
s'écrit = 
NÉE ES poule ae > 
lc pre À [2moc + U+ Zi] : ÆS 


X [(p+U}+/%(RotU.c)] 


_ ékI(grad(Wo+ H).e][(p+U).9] . 
[2m0e = U6= if les — 0. (18) 
Posant S 

I =p+U, (19) 
cette équation devient 
Fee ee 
le re 
Emo CUS Eh (RotU.c) 
AD 


[emuwc EU +kP 


x | (grad(o,+).n) 


— ife(ead(un)nn)]]}es= 0 (20) 


à Eù tour si l'on à M = rot U re tondaie 

par exemple au cas électromagnétique U = €, 

— 3%) on retrouve en première “Hauts 

set supposant le couplage électrique faible, U,< m,c, 
14 équation de Pauli 


re ; [p + RAP (ae a) lei = 0. 


l Ci È 2 m0 (2 20 


Dans le cas des interactions centrales, U rentre 
dans le terme de phase, et il nous reste à partir de (20), 
introduisant le moment orbitral 


7 L=x/p, 
LE--men—ntr) 
—. p° 
PSE ve ne Ti] 5 
TENTE FL, P Es PS ape 2| @s=0. (2 ) 


Nous obtenons ici pour % une équation qui 
_met en évidence les termes de couplage spin-orbite 
bien étudiés dans la théorie de l’électron (1). 

Ces termes n’interviennent pas dans l'hypothèse 
de faible couplage, auquel cas nous sommes ramenés 
à l’ équation de Schrôdinger 


E ; 
rec P Je (22) 


20 C | 


_ Un cas particulier remarquable correspond à 
. la compensation exacte de U, et de J, 


a Le à 


no ce cas le terme de couplage spin-orbite est 
écarté et il nous reste l'équation de Schrôdinger 


re p° ee (23) 


Hi (Soso+(S.5)) + Las. 


Nous ne considérons d’abord que le cas de l’interac- 
Étionen Las. 
Les équations en 9,, % 


Foret 
© 
Ÿ 
(= | 
Q 
: 
ie 
| 
se 
TK 
PX 
Le) 
A 
nn 
À 
LI 
NS 
9: 
KB 
S 
(l 
© 


nous donnent immédiatement sous la seule hypo- 
thèse E < m,c?, 


(grad. s) fe = (25) 


Btc 2 Mo C MVC 200 


(4) Voir, par exemple, pour une discussion de ces termes, 
- CoNpoN et SHORTLEY, Atomic Spectra, 


L en ae cuiée Sels (r), 


CHERE ME TION 


PC 2 Mo C 2M0C 


(grad,./:) (==) Fes = 6! 220) 


20 C 


Nous sommes ramenés à une équation de Schrü- 
dinger avec deux termes d'énergie potentielle, le 
premier 1;(r) n’interviendra que si l’on est dans un 
cas de fort couplage et ce terme sera seul à inter- 
venir si nous avons /, — const. Le second 


7 (grad, 1)(ES) 


2 Mo C 


dépend du spin et introduit un coefficient d’interac- 


tion analogue au magnéton de Bohr [on peut écrire 
1 =9g#,(), #3 étant le champ pseudoscalaire 


extérieur agissant sur le corpuscule caractérisé par 


his ; 
la pseudocharge g. Tr Sera un pseudomagnéton 
20 


relatif à l'interaction pseudoscalaire. 
Si nous considérons simultanément les interactions 
en S% S et Z,, nous avons les équations en 9,, 


E va 
om HAE (S.c)|: 


Hp ES belge, | 

E | (27) 
Ê ad) ) es + [(p.0) + Si— 12 ]@1= 0. 

Si nous supposons E Z mc et le couplage en S faible, 


S|m,c, nous obtenons pour w,, en introduisant 
le moment généralisé 


= p.506, (28) 
CELA STE 2 
le 2 Mo C 2 M0 C 


ae L(s SE 7 grad/7; 1 | lee 0: (29) 
20 C 


L'’interaction en S, nous oblige à modifier la défi- 
nition de l'impulsion, d’ou l’expression classique de 
l'énergie cinétique, mais n’introduit pas de terme 
d’énergie potentielle dépendant du spin. 

Si S, n'intervient pas, nous avons deux termes 
d'interaction en r. Si en particulier S est le gradient 
d’un champ pseudoscalaire (cas du champ mésique 
pseudoscalaire) tel que 


S=/fS, B= 8% avec $ = lo gradÜo, 
nous obtiendrons l’équation 
(EEE RE 15 


le 2 MC 2 M0 C 


+- (Fi À 5) (gra os) les — CN (el) 


270 c: 


Si en particulier #, — %,(r) le terme d'interaction 
< \ X" 
se ramènera à la forme grad,.#, (==) à 


On voit nettement sur l’équation (29) l’équiva- 
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lence déjà signalée par plusieurs auteurs entre les 
champs en Za; et (S. co) à l’approximation non 
relativiste. Toutefois cette équivalence n'existe 
que si 1? et S sont négligeables devant m,c, c'est- 
à-dire dans le cas du faible couplage. ë 


30 H=(P.r)+(M.p). 


Nous avons déjà discuté dans un Mémoire anté- 
rieur les interactions de ce type dans le cas où 


(P,M)=/(6,4€), 


&, #€ étant les champs électromagnétiques. 
Les équations en 9, et 4, s’écrivent 


ET ce = +] c1+[(p.5) —1(P.0)]2: = 


E % 
[£ = (M) 22 


Si nous supposons E < mc, |M|—<m,c (faible 
couplage magnétique), Cnrogaction de l'impulsion 
généralisée 


(31) 
+ [(p.7) +i(P.0)]%1= 0. 


= piton) (32) 


nous donne 


fe M =) ji pr+%divP|ln— 
Le 2 10 C 


(38) 


a considération d’une interaction en P introduit 
nécessairement une impulsion généralisée IT dépen- 
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je Ne 40. 

4° Si nous examinons les uote ( 1) + A2. 
nous voyons que If ne dépend du spin qe par 
l'intermédiaire de P et.de Sy. 

Si nous écartons les difficultés relatives à cette 
dépendance en prenant P — S; —o et si nous 
supposons toutes les interactions radiales 


Uo = Ui(r), LT B= (Tr); 
S 35 
S=<S(r), M= > M(r), U => Ur} Ca 
l'équation (12) s'écrit , 
FAURE he # 
- . 2m0c ARR 2 M0 € PA 
+[sO=MO EEE À = grad 1e 5 Et} 0 CR 
ou encore, posant ; 
15(r 
TO TO 1 = vtr) | à 
(37) 
Sr) M(r) F2 — grad, +0 = F(r);. 
FA Apr e en ee FEAT 
Le: one. For) | F(r) A à Es res 


* Passant en coordonnées polaires avec 


= rsin0 cos», y =rsin6 sino, 37 COS HS 


cette équation nous donne Aout les fonctions Pa (pans : 
%29), le système 


dant du spin. Par suite, on doit en conclure que dans E TRS A LR 
le cas particulier de l’électron possédant le moment Cr oi | Tr a | 
magnétique intrinsèque de Schwinger, il est néces- RSS EE RS TE 0 
saire de modifier la définition de l'impulsion, notam- [rome r° Ée sCoet 98 
ment dans les calculs de diffusion coulombienne. 20 dou 
Il en serait de même dans l’étude de la diffusion cou- | RE | AC ae RE g 
lombienne d’un nucléon possédant un moment eig sin8 vos + COS 0 voi = Ë 
magnétique intrinsèque. AN nt C2 
’ X * 
Dans le cas coulombien P —P(r) ms PU) =) Posant 
l’équation (33) se réduit à o1(7, 0, o)=eimpon(r, 6) 
ee (M.5) gas = ein tlpess (7, 0) (40) 
Le .(m entier), 
GS. è 
o Q == Dior 4 É x 2 Ê 
270 C pe se )( r )| He système s'écrit encore 
OU ER ae nr Ne nu | 
: EG >M0C | Or?  r or Fur Don 
LES CA RCE 0 me? - Ë : 
RE ne SU = + ?, S Do ES S ( Door À 
Es r? É soie 5 5 | ARS Rene 0 Ces 
—+ 
R2: NT; T0? d (m—1} = l DU 
—— - — Gh ? do + F S (001 —= 06. 
| Re É cotg s0% RES | 3) (0080 | #22 3 (r) Sin Ü 291 = 0 
III. — Nous allons maintenant examiner plus valeurs de trois intégrales premières : l'énergie W, 


en détail les conditions dans lesquelles l'extension 
des méthodes de Dirac (1) et de Darwin (#) permet 
la résolution de l’équation de Dirac dans le cas 
d’une interaction générale. 

La méthode introduite par Dirac repose sur la 
caractérisation des fonctions propres 4 au moyen des 


le moment cinétique total 


(49) 


li l 
NEL+E Se (p)ESS 
et une intégrale particulière remarquable \ 
u = a [(Li3) + #7) (43) 


n HO MR, Late 
x L 
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que nous ne leon intégrale première de Dirac, 
qui permet notamment dans le cas de l'atome 
d'hydrogène de séparer les fonctions d'ondes en 
produits de fonctions radiales et de fonctions sphé- 
piques. 

Nous allons donc d’abord examiner ici dans quels 
cas nous avons encore ces intégrales premières. 

- Nous considérerons d’abord le cas du moment 
« cinétique total N. 
On sait que l'équation d’ ondes étant écrite sous 
la forme e 


+ 


s : RU x) + mo Cas + Hi]Ÿ = 0, (44) 


pour que là grandeur correspondant à î opérateur C : d'ondes on remarque que l’on peut écrire 


soit intégrale première, c'est-à-dire telle que 


4 Cu dd 210 
a) * pe 


on doit avoir entre les opérateurs Cet £ la relation 


[C, ÉI=C£=£C=0. (45) 


On voit immédiatement que 
ki He 
L - 3 % Po (p:x) + MoCas 10 


Nous sommes donc amenés à chercher dans quelles 
conditions nous avons | 
= [N, io e [Do + (UÜ:x)] 

+ ((P.x) En (M:w)) + (So (S.s)) —- L2;] 10: 


__ Considérons un des termes de Æ,, soit Ax,, nous 
avons 


en 
= ih[x/\grad A]as+ ŸA[oas— aus]. (46) 


_ On voit immédiatement que cette expression sera 
nulle si nous avons . 


HA Po) Pi 02 F3 CE x/\grad A = 0, 


condition toujours satisfaite si À = À (r). 
. Plus généralement, soit À un vecteur radial 
MA OT. 5025200) 


| PE 
(A:x1) =(A.6)ex= A(r) (E jee 


Or, on voit immédiatement que 


à CRE … A(x) ee 
er [(x/p ); (As)pe] à cl A ook, 

le A s 

Île, (As)pr] = = © ie (x.c)l64 


DAC) 
= —i% D) ixA or 


d’où 
2 + Fa (assa]=e (41) 
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Par suite, le moment cinétique total N sera une 


intégrale première dans le cas où A, aura la forme 
générale 


Hi= E(r)e + Uol(r) + Ur) (Æ FE) 


POS) M (EE) 


+ So(r)50 + St (£ =) + Rr)e. (48) 


Nous allons maintenant chercher les conditions 
relatives à l’hamiltonien d'interaction pour l’exis- 
tence de l'intégrale première de Dirac (43). 

Pour mettre celle-ci en évidence dans l'équation 


EE» EE) 


= (ap + (2 (==). «) 


7 


L’équation (44) s'écrit alors 
Pre) 
i(Æs) (EE) © mens ]y= 0 (60) 
et en introduisant l’expression (43) de u | 


BOIRE 


(=): NO CA Hi |4 = 0 ÿ (51) 


Fr 


Pour discuter les conditions dans lesquelles u est 
intégrale première, nous calculerons le commutateur 


[uf — Lu]. 
A partir des expressions 


æ[(s.L)(p.x) -(p:x)(o.L)|+2/%x;(p.x) 
= p:p1l(s.L)(p.5)+(p.c)(e.L)+2#(p.0)], 
_(s:L)(p:5) +(p:s)(s.L) =— 2h(p.«), 


on voit immédiatement que 


[u, Po (p:x) + moca]=0, (92) 


u est donc bien intégrale pipes en l'absence 
d'interaction. S 
Pour discuter les conditions de commutation 
de u et de H, nous distinguerons deux cas : 
a. Fe HE A ph 
{ pspexl(s.L)A — A(o.L)] 
CRE O2), | 
pelCs:L)A + A(s.L)+2#%1] ( 
È (Æ=1, 2). 


Prairie (53) 


Remarquant que si À = À (r) 


LA(r)— A(r)L = 0, 
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on voit que u sera intégrale première dans les deux 
cas 


Hi Ui(r), Hi =) re (54) 
b. Hi= (ASv)er 
Nous avons alors - 
ospel(s.L)(A.s)—(A.c)(s.L)] 
(KES 0 0): 
us Hu sox[(o.L)(A.0) (59) 
+ (A:6)(s.L)+2/(A.c)] 
CA; 2e / 
# 
Or, on a 
…(o.L)(5o.A)+(s.A)(o.L)--2/(6.41) 
= 2{[(A/x).p]+1.(x.RotA) 
+ À (grad(x.A)— x divA +A).c|. (36) 
Par suite, si ss 
FE se 
A—= à A (nr), d’où Rot A = 0, x/\A = 0, 
pour À = 1, 2, on aura 
uH,— Hu = hosex| (grad(x.A) 
| — xdivA-A).c|. (57) 


Tenant compte de 
(x.A)=rA(r), 


grad(x.A )— 


div A = 
| A(r) 

XX nt 
= p 


+ 70) 
nous voyons que u est intégrale première pour 
I= AG) (ES); =). (58) 


SIkR = 0,67 ON à 


Hi= P(r) (x 


Lu, Hi]= cer [én(x Rot A) 
+i(o.[2(A/L) 
—ili(x grad À —x div A)])|; (59) 


mais (A AL) ne s’annule pas en général et notam- 


ment si À — À (r)E, il ne peut y avoir compensation 


dans les termes en o et par suite dans ce cas [u, H,] 
ne peut s’annuler. 

Nous avons donc mis en évidence a cas dans 
lesquels l'intégrale première de Dirac existe. 


19 IT, == Uotr), 

C'est le cas considéré par Dirac (atome d’ hydro- 
gène). 

920 Hi = trier. 


C'est le cas étudié par Caldirola. 


3o H=P (EE). 
Fr 
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‘ li 2 112 
"u?0 = [(- 59) — | | 
' 9 4 


ce qui montre que j ne peut prendre la valeur zéro. 


C'est le cas considéré par Breit. 


Ce cas trivial peut être écarté par modification 
du facteur de phase de la fonction d'ondes. 
Pour utiliser cette intégrale première, Dirac a 


montré que ses valeurs propres sont tous les nombres 


entiers sauf zéro. 
Le facteur «, ayant pour valeurs propres Æ 1 et. 
commutant avec (L. 9) +; les valeurs propréss 


de u seront + celles de (L. o) + . D 1 
ub = a; [(c. LOS TRUE ; Fe 
Tenant Copie de l'identité 2 2 : | 
(s.L)(s. Dane 25 | LS a 

s à 1 o? 

= L?—-—"/>? 5 - 7 cot£g 6 - CU RE RE de? D can) à 


posant SES 
Ù a Pr #) eirro Se . 


ES 


en. 
RULS 


Par suite, les valeurs propres de É D: ne. F et. 
— (1 +3:)# et l’on a soit l ne S-> 


J=Æ4, 


où P/' est le polynome de es associé, on 


déduit 
(s. L)[(s. Lo 4 = —= RO 


Soit. 7 = Sr + 5 


Mais, d'autre part, on a 


= a[(s.L)+%] = «| [<.(L 2 


Donc j prend toutes les valeurs entières = o où <o0 o 
sauf zéro. À 
Dirac à remarqué que les trois matrices (& 0). 
(x.ioa), «, forment un système de trois matrices . 

unitaires anticommutantes, car l’on à 


X.a \? 
La IR 
F 


(=?) ee in Se RS 


Par suite, si la solution de l'équation d'ondes 69. 
dans le cas ou u est intégrale première, est telle que 


uY=.jY; (67) 
on pourra écrire 
d=y(r)e8(6, 0) 

et posant = 
Ne ; X.laœu; ; RS 
(É)-4 CES 
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2 


Hi= Ar) Sa  : 


s 


y 


bé 
> à 


X.1a; \? 5 TMS 
—— | —=1; - LS 
r ; € F'ASES 
X.x XX. 9 
RS r Ar Or, Creee ne =—=10: EH 
r v ; Fe 


D'autre part, ces matrices commutent avec © 


207) sû sera solution de l'équation 


[w 


Five 
d se pr Ft ei 2 + Mo Ces 


+ H(r)es Da(r) + P( r)enVx(r)=0 (69) 


& par suite, les Fctions d'ondes 4, ne dépendront 
que de deux fonctions radiales 


à 


xa(r) et 207). 


- IV. — Nous allons maintenant examiner dans 
uelles conditions l'équation (69) admet des solu- 
ions dans le cas généralisant celui du en 
oulombien, où nous avons 


Los dis Do ge ce a ant des constantes. 
4 Nous utiliserons pour cela une extension de la 


(71) 


Moc = Ho, 


r 


épendant des matrices p/, p,. Pour r > cette 
quation se réduit à 


He 0) e : É 
fées ge + eat dico piGn+ bn] #07) = 0, (73) 


4. 2 cù — (ue Go) | +7) = 0 (74) 
ÉRons poserons 
(75) 


a? ae ue bo} + cÿ — (æ —- ao Ÿ ; 


6? > 0 correspondra aux états liés de la particule, 
2B2<o aux états libres. 


_Si«8? > 0, l'équation (74) écrite 


2 
É A af | e(r)= 0; 


admet la solution tendant vers zéro pour r -> 


Ç 


o(r)= À e%8r. (77) 


(5) TEMPLE, Proceed. Roy. Soc., 1930, A 127, p. 349. 


s (0 . 2 es Jet) = 0 cm) 


(r) représente ici un couple de fonction 9,0%. 


à laquelle correspond l'équation du second ordre 


(76). 


“SUR LES ÉQUATION | D *ONDES. DES. À CORPUSCULES DE SPIN #/4x , 21 


Nous chercherons alors une solution de (72) en 
remplaçant À par un polynome en r, soit 


2° rs e—ahr, (78) 
les constantes À, sont ici des matrices | se , 


Introduisant cette expression de  (r) dans (72) 
et prenant le coefficient de r°-!, on obtient immé- 
diatement la relation entre coefficients A. 


Éd tpes — pi — a) + bies]As 
+ [(w + a6) — 10, 28 + pr co ps (Hot bo)] As 1 = 0. (79) 
Si pour un indice s, nous avons À, —0, 
avec À;, ;< 0, nous devons avoir 


[ai tp250 — pi — 1) + bip3] As, = 0; (80) 
ce qui exige la relation 
É Lai LS (j— )}—0bi]4,=0 
et, par suite, 
Ft ee (81) 


Multipliant (79) par 
= 1pa af — pi co — ps( Mo + Vo) 


+ Ag 


et tenant compte de la définition de «6°, on obtient 
entre À, la relation 


[+ ao+ dpa ap ps (o + bo)] 
X [ai doss — pit — 1) + bips]As— 0. 


£1Co 


(82) 


Si la série des À,, s’arrête pour un indice s’ tel 
que A: — 0, on aura pour cet indice 


Le + ao tp af + paco + pa(to+ bo)] As = 0. (83) 
Tenant compte de cette relation dans (82), celle- 


ci devient 
ai(w + do) — s'aB +(j — ci)co— bi(uo+ bo] A =0 (84) 


et, par suite, s’ sera déterminé par la relation 


s'ai = a(w À &) t (7 C1)Co bi(o i bo), (85) 


qui s’écrit encore 


[ai (æ + ao) + (J — Ci)Co — (10 + 
(Bo bo) + ci —(w + mm} 


. bo) b1 P 
“1e : 


(86) 


Nous déduirons les valeurs de W de cette expres- 
sion en écrivant 
S = $ +2, (87) 


n étant un nombre entier et s, étant déterminé par 
la relation (81). 


On obtieut ainsi l'expression générale des 
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niveaux W par 


a [(uo + 0o6)b1— (7 — ca) 


ee : Fee ai (son) 
I 
- ai +(ss Fran}? 
< nn + n}? [ai (« uo + Do)? + cè) 
Re —((j— ci) cote bo) 1) ] 
c = \ 1 
+ (son) (( Mo + bo)? + ci)] fe (88) 
LHC A partir de cette formule générale on obtient 
facilement les niveaux dans des cas particuliers 
intéressants : 
10 Cas électrostatique : U, seul, 
aÿ(æ + &o)° à 
2 — EE = 12 —= 47 
(So+ nr)? = Be a , SO = ai 
# (55 > 0), 
ch É 4 
as r D + do = Ho 1+ Se TA ANR S (89) 
| (n = Va?) 
FR < qui pour 4, — 0, nous redonne la formule de Dirac. 
= À l’approximation non relativiste cette formule 
: : nous donne 
2 ‘ F ie E =w— mic? = — Le: EE (90) 
| 2 n +) 
2 bi 
2 Tr) = | b, + Fe seul, 
Re 27 b? (to + bo }? » à o ) 
* (Son) = TA TR EN E SG =J*+ bi, 
ï 1 
SE ; œ = (Lo + O5) Ç (91) 
mean 
: “et en première ROpORNNaton, Si. —0, b,ipetit 
ju devant }j, 
AS JE re |" ; h: (99) 
2 En + | 
expression analogue à celle du cas coulombien. 
Nous retrouvons ici les formules de Caldirola. 
30 Si nous avons simultanément les champs en U, 
et, J, avec 4ç — di ==.0 et ai = bi, d'où s$ =|j, 
posant n +|j| = N on obtient 
à M a; 
Si &i= bi, W = — Re 9 
: ai 1; mo C N + a? (93) 
SD 1&1 —1 = 01, W= = moc? USPRN OR (94) 
‘ F IN? == Ph : 
A9 P(r) = [co 1x 2] seul, 
ue (GESXBECT pire 
PR le où nu Wed Ne VAN 
d’où 
Prius 6; es PR (95) 


n+ly—oc: 
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. 4 sE k LERS Re LD 40. 

| | Lie “A 
Contrairement aux ca précédents, l’ existence des 

niveaux est liée ici au terme en (Ar 
Si € — 0, nous obtenons simplement 


il 


lice] c 


os nie (2e )] (06) 


et à l’approximation non relativiste 04 
h2c$ ( Î j 
2mo \n—+|J| 


seul ne contribue qu’à un dépla- 


Fe (97) 


* CG 
Le terme en : 


cement des niveaux PRE dans les cas pré= 
cédents si 4, “0 où b,-<o en remplaçant alors BE 
par j — Ci. 

C’est ainsi par exemple que, dans le cas des inter: 
actions électromagnétiques, nous obtenons si une. 


interaction en ss due à un moment électrique 


: 15: : ie (22 
s'ajoute à l'interaction coulombienne en Æ 


ai | £ 


W = moc°| 1 + RTE DR) 
= Ga a) 


Toutefois nous devons remarquer que lines 
duction du terme en c, lève une dégénérescence sur ne 
niveaux. Ë 

En effet, dans le cas électrostatique coulombien 4 
on avait 


S 
L * 


(98) 


Hahn cw 7 


j n'intervenant dans la formule (89) que par son. 
carré. 
Au contraire, si le terme en E intervient, nous 


avons 6 : +, F8 
2e 
W; =; 10 C2 L ae ( ES | % È 
n SE V( Ca) s 
VCJ t) net ae (100) 
We; =—"}n0 C2 |: NS APR en eee ne £a | PE 2) 
(nt +a)— ai) } 
et en première approximation FER 
Pi W; Se 2C1 a} : (100) 
/ 6 . IQ C? (n + PARL 


Ces résultats peuvent être utilisés dans l'étude des 
niveaux d'énergie de l’électron dans le champ - 
coulombien, lorsque l’on tient compte de l’existence 
du moment magnétique intrinsèque introduit par 
Schwinger, d’où le terme 


Le | 
HI = (6;7), 
ACTES 
& élant le champ électrique correspondant au 
potentiel coulombien 
&(r)=— 


No] 


) 


à A PA à « J Es SA c 74 
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: d'où Rec £ de la représentation des matrices « développée par 
2 CRUE mer (aen). (102) Franz (16) et par Sauter (1?). 
4 | À cr? Le : Avec la représentation des x considérée l’ équation 


. Mais ce terme introduit dans l’équation d’ondes 


MT . CP: I 
- radiale une SeMrRe en — 


NW On peut éviter celle-ci en remarquant que dans 


- le cas électronique il s’agit d’une perturbation très 


Re et que, par suite, on ee en PEAREE approxi- 
. mation remplacer le terme en …. ; par = —, où r, est le 


_ rayon de l’orbite de Bohr noce au en ibie Les 
pique N =n+ij|soit 
FRA LE 

= 


/ moe? 


 Posant p. — pu, où À est une constante, 


Re eh SPACE 
4 UT moe ER TE 
on voit que 
7: 5 ka? 72 
a ] Ci =— RE Ar AZ, 
: d’où PC 
ee W;= W ; xai Zi 
- A; = ASS E PRET EE Se £ 
7 Mo C2 (n +171} 4 
Selon J. Schwinger (f) on doit avoir À — — > > d’où 
re 5 | CAT 
De. A — 104 
EE - DRE CES) (10) 


_  V. — Nous allons maintenant examiner de plus 
= _ près la forme des solutions de l’équation d’ondes (1) 
dans le cas des interactions radiales en U,(r), L(r) 
- P(r) en utilisant la méthode de Darwin mais avec 
la représentation utilisée par M. L. de Broglie et 
introduite en (1), (2), (3) pour les matrices &. Il 
. serait toutefois possible d’exposer les résultats qui 
suivent en utilisant le formalisme indépendant 


. ar : f of ke u+1 
2 rai) Aer PT Ce ne 
à Quid) f(r) FE = . (Lin re 
: Q (22 
HITS — dr 


r 


Suivant la méthode utilisée dans le cas de l’atome 
. d'hydrogène, nous cherchons une première solution 
de la forme 


7. W= 104 F(r)FE,:; 
14 F(T)VE 1 


k +17 


UE = d; 


GRIS) 


‘10S 
di = dr, G(r) Fe { ( 9) 


-  Reportant ces fonctions dans l’équation (106) et 
utilisant les relations de récurrence entre poly- 
nomes P/(0), on est conduit à un premier système 


d'ondes qui s'écrit 


[TE + Uo(r) + (re —— Li(r))e 


C 


5 + (pioer = PG(EZ ) 2 [= 0, (105) 


se développe suivant le système 


[E À Us + Mo € + n | | 
C 


+ih[(dr+10;)b;+d:43] 


EI + )h+ asbl 20, 


W 
[= + Us + moe + 1 | Vo 

+ 1h[(dx—i07)b3— 0,4%] 
zh] 


= 02 


EE pa 
Se (106) 
É + Uy— (moc + Hi )| UE 


+ ERT(dx + é d;) Va +0: di] 


10) 


Suivant la méthode de Dar nous cherchons à 
représenter chaque d comme produit d’une fonction 
radiale et d’une seule fonction sphérique du type 


+ Ia + iph+ au] 
[= + Uo— (moe + I )]# 
PA = 10;)di— db] 

+iË te y) hi] 


CRU L £ d k+u s 

M ———— eiu? sin 0 | —— D -EIX 
y} Tee ® sin (rx) [cos20 — 1] 

— HEC ue (07) 
Nous utiliserons les relations 

Guen (Es Ep) rest 

KES \se 

À (RE &)(E 4 ù 0 (LE - = y) ne (408) 


7 


d df . _k+ 
pus ) VE sd (x + u)( ka) (+ : EE 


de solution qui s’écrit, à un facteur de normalisation 
près; 
—1FYY 


k+1)2 


D=(f— ur) G 
Ve PE 


d, = (Æ 


D 
= u)G x Y to : ) 
(45) FrANZ, Sitzungberichte München Akad. , 1995, p. 380-435. 


(7) SAUTER, Zeïtsch.f. Physik, 
P. 299. 


1930, 63, p. 803; 1930, 64, 


1e 


F (r) et G (r) étant solutions du système 


PEUT 7) me + AD] E NEREN 
. DR LR PT) RU 
4 ; ose ñ Je = 
> 1 [W 
; x [= Ur) (me nn)]e 
NE ÉTÉ PO] ED 
ÿ de TAC O TES ° 


\ 


De même on nee facilement un second système 
de solutions | 


a. 1 W= if u)FTE ss, 
Dre tre ES u SEE, 


Ya a GYF, 
h=— GE, 


LE Pa @) et G(r) étant cette fois solutions du système 


PL 


Ar re. ; 
SRE RE Ro Ad Pet 
ne ë ARE Are st = ]e=c 


les ee (110) CEAURLE) s’écrivent encore en 
introduisant la fonction v;(r), telle Èe 


; pi=f; Pr "8" | 
À LOW + Da) + (moe + H)es] EE 
et) ARE ET NC Re 

A rm 


mA 


7 IO + U)+(moc + Tips] 


“ & P 


Fiorse Sr pjeimee. (113) 


le système (69) sous la forme (72), avec pour On, 
k=ÿ, pour o5 k +i=—7. 


Nous allons examiner la forme explicite des 


solutions de ce système dans le cas 
à 5 Ci 
Us = const., li = const., P=——. (114) 
# 
UT = , Nous obtenons alors immédiatement par itération 


Le NAS 
LIL RUE TS ANT 


U) (Mo ce + | 
| 


=0 (115) 


r? 


compte des équations du pret miel 


TA Ne RU uours - 


Nous retrouvons donc par la méthode de Darwin 


SR 


_cune de constantes à qui ne 
dantes et pour les déterminer 


A partir de ‘ De nous écrivons > 


£ { | k Fe Ÿ ï F = Fe u . 
— +Uo+ mc + | 
tas Re 


TR 


Le hote 


LS 


et utilisant les formules den récurrence 


é = 
ar" &n) LS Ka ? 2 È 


Rs °. 


AA ire Z 1) 
Rues 


on trouve facilement 


An —= K.. 


1 On en déduit la Non on 
one 


+ BR CE) 


vr 


ne 


k—1) 


PPT ANT à 


2 


7 :(K 


ÿr k+e+ 


TR . 


A 9 
Cr 
> 


DEL 


LS: A) 


AE A 


12 


192 
ETAT ( ) 


De AU 
+B TZ L— 
21 PE 2 . 


Vr fers 


Kr 16 
Pret ) k 


TE 
1 (Kr)YE u—1 


Pb,  (Kr) FE. 


=Z HAS 


> On Gbtient ainsi une solution générale qui dépend 
| de quatre constantes arbitraires. 


VE. __ Nous Allons maintenant revenir à notre 
équation générale (4) et examiner si, en dehors des 
cas dans lesquels existe l'intégrale première de Dirac, 
il ne nous est pas possible d’obtenir des solutions se 
… présentant comme des produits de fonctions radiales 
. et de fonctions sphériques. 
Le cas le plus simple dans lequel nous nous pla- 
_ cerons ici est celui des interactions d'intensité 
_ constantes. Nous allons montrer que dens ce cas 
il existe pour certains types d'interaction une inté- 
. grale première qui nous permet laséparation cherchée. 
En effet, dans le cas d’absence de champ extérieur, 
l'équation de Dirac 


[Po+ (pa) + mica]4 = 0, (193). 
à admet notamment l'intégrale première (18) 
ï. v —(p.s). (124) 


Dans le cas où nous complétons l'équation (122) 
par le terme H,, avec des termes d’interaction cons- 


_ tants, on voit facilement que v est encore intégrale- 


_ première si nous avons 


"# Hi= Tia + Uo- So €: t Lp, (125) 
@ avec 


Le Us, So; 1—= const. 


Dans ce cas nous avons également encore l’inté- 
grale première du moment orbital 


li 
N=L+ 6. (126) 


(5) ProcA, Thèse, Paris 1932; G. PETIAU, C. R. Ac. Sci., 


Paris, 1946, 223, p. 269-70. 
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Soit 1% une valeur propre de v et 4 la solution 
correspondante, telle que 
(p:s)Y = AY. (127) 


L’équation du second ordre que nous dérivons 
de (123), complété par (125), nous donne, tenant 
compte de (127), 


(128) 


D'autre part, l'équation d’ondes se décompose 
en deux groupes de deux fonctions ol}, 0°), tels que 


FE E Uo+ moc 1] (1) 
+ [(p.s) + So+ih]el 
[+= (met AI 
A [p.c) + So—irijot= 
Tenant compte d2 (127), nous voyons alors que 


les deux fonctions ol) seront déterminées à partir 
des ol) par : 


= O, 


(129) 


À li + 90 EE AE 


se + Us + mo + Ji 


Il nous suffira done de déterminer les solutions 

de (127) avec les valeurs (128) de 1%. 
Or, les solutions de (127) satisfont à l’équation 
(ps? Ÿ = pv = 22%20. (131) 


Soit y(r, 0, ©) une solution de cette équation 
qui s'écrit encore 


d? HATU Es 
dr? ror 


1/08 0 : 


Posant 


I ms) Ce 
sin2lo do? KE 


= (Cr) etimp Py(0) sr" 1(r) F}(0), 


on voit immédiatement que 7(r) sera solution de 


9? 2410 ACER) * > 
É  r or r? rx ut) 


Gr), 


(133) 
d’où 


RS 
{(r)= Ar ?Z (134) 


1 
l+; 


A étant une constante et Z, 10) une fonction de 


Besse] d’ordre demi entier ! + : 


Les fonctions 4°) se déterminerons à partir d’une 


solution y*”\ par 
CDHALAEE [(p.s) + Fe R]xbrm À, (135) 


[L + | APTE — Ce 2 Jens, (136) 


(360 nous donne immédiatement 


— NT. ex Ge AG FE “ 


pig De DS 


Tenant compte des relations de ni or. 


PR: 


: : de. ] ne IP Er F, pret 


Or ES me 
a 7 Le 


Guide, 1 nr 


à Mie , Er Eat E CT 

à 5 ES S < É $ L ; l'a PES #16 Ÿÿ FRERES + 
dan, none nn Le 
PEL RAR Li ai s FT NU A PET 
Nous en déduisons la solution A et RE ne 
30 ? L #0 LD - PE Dee # LA: ï 


F 2 TT 


ge ha = ES 


FLE 


CE RES 


Lee 40 r). FE m2 m 0 


se an Er = CUITE mn 


RE REEEE 


san, | Uementen, 


Le 
JS 
2) 


k,300) Sn ns 


ue 


e “ À i £ . ne à 


Age ae br de À correspondent eux 
_ solutions de ce type pour les fonctions œl1 , les pen 
Ps re en déduisant par la relation (130) et l'on voit. travail = 
Je encore que la solution générale SUN de share RIRE 
Hours arbitraires. RONDE NES Te md es DU RC NPA CAE 
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SÉRIE VII, TOME X, OCTOBRE 1949. 


ÉTUDE THÉORIQUE DES IMPULSIONS DE TENSION ISSUES DES CHAMBRES D'IONISATION 


Par P. BILLAUD. 
Laboratoire de Physique des rayons X. 


Sommaire. — Au moyen de fonctions représentant le mouvement des ions dans le champ il est 
possible d'étudier théoriquement les impulsions fournies par des chambres de géométrie quelconque. 
Certaines lois sont établies, qui permettent de tirer de l'examen d’une impulsion des renseignements 
sur la particule ionisante, notamment son énergie et sa direction par rapport au champ. 


TJ. —— Calcul de l'impulsion. 


À. Mouvement des ions. — On sait qu’une 


. particule chargée en mouvement crée dans un gaz 


_ des paires d'ions; ces ions se meuvent dans le gaz 
sous l'influence d’un champ électrique en cédant 


constamment l'énergie que leur communique ce 


champ aux molécules du gaz, par chocs. Nous 
 retiendrons surtout le fait que cet échange d'énergie 
est, en moyenne, irréversible, et que le travail 
… accompli par les ions au sein du mélange gazeux 


F 


. constitue toujours une perte définitive pour l’énergie 


a potentielle totale du circuit électrique relié aux 
. électrodes entre lesquelles était établi le champ. 


"B. Équation du circuit théorique. — La 


chambre d’ionisation ayant une forme quelconque, 
_ nous supposerons que son électrode dite « collec- 


_trice » est reliée directement à la grille de commande 
d’un tube à vide de type quelconque et que cette 
grille n’échange aucune charge électrique avec les 
autres électrodes du tube. Le conducteur composé 
de l’électrode collectrice et de la grille est relié à 


- Ja masse (sol) par une résistance pure À. On applique 


à la chambre une tension V, positive ou négative; 
ses deux électrodes présentent une capacité 
mutuelle C,. Pour calculer la différence de potentiel u 


_ à laquelle se situe l’électrode collectrice par rapport 


au sol, lorsqu'une ionisation transitoire est créée 

dans la chambre, il est nécessaire de tenir un compte 

exact de toutes les influences électrostatiques 
auxquelles elle est soumise. 


On a schématisé sur la figure : : 


La protection électrostatique de l’électrode collec- 
trice (blindage à la masse) et la cathode de la lampe 
d’entrée, par la capacité C4. 

La grille écran de la lampe d’entrée (s'il s’agit 

d’une tétrode ou penthode) : capacité C;; potentiel 
fixé, Ve. 


La plaque de la lampe 
variable V, = x — fu. 


: capacité C», potentiel 


Un conducteur quelconque susceptible de pré- 
senter une capacité mutuelle appréciable C;, au 
potentiel fixe V.. 


7 


diverses dérivations des 


les 
los las etC., orientés comme sur la 
figure (1). Les capacités sont chargées à des tensions 
instantanées 


Il circule dans 
courants 1, 


Fr 


Vo = Vo — u, Ve = Vr— U, Ve Ve— U, 


Vy = — U, Vp=a—(B+1)u. 


On a, d’autre part, 


(42 


R 


= = 19 + lotte + 1 + pe 
Soient w(t) le travail effectué par le champ de la 
chambre sur les ions à l’instant t et p (0 — su la 
ac 
puissance correspondante; on peut écrire que la 
puissance développée par l’ensemble des forces élec- 


tromotrices est égale à p ({) augmentée de l’accrois- 


(:) Le courant ( est un courant ‘partiel qui se superpose 
au courant de plaque proprement dit, de même à. 
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sement des énergies emmagasinées dans les capacités 
et de l’effet Joule dans la résistance R 


Vo io + Voie + Vete+(a — Bu)ip 
= p(t}=— Co(Vo—u)u — Cx(Vr—u)u 
— Ce(Ve— u)u + Cyuu’ 


du 


L ; u? 
—(B+n)Cpla—(f+i)ule + avec Ur ae 


Sauf i, et à, les courants sont les dérivées des charges 
des capacités et l’on a 


SOU UNE AC EUR —=— Cyu, 
ip=—($ +1)Cpu'; 
d’où 
6 u \ n S 
= = +[Cx+ Co+ Ci+(B+1)Crlu 


R 
En portant ces valeurs dans l'équation précé- 
dente, et après réduction, il vient 
(Po— u) +(Po—u) 
SR [OC + CEE (BE DOru = p{i) 
que l’on écrit, en posant 
CC DCE OT be D Ce: 


CIO 0 


Dans de nombreux cas expérimentaux u est 
négligeable devant la tension V, de la chambre, 
et la relation ci-dessus se réduit alors à 


(bis) 


Nous nous limiterons à l'étude de cette équation 
approchée. 


Remarque. — Lorsque la résistance R est très 
faible, l’approximation ci-dessus envisagée est valable 
a fortiori et l’on pourrait étudier le courant t, ou 
mieux le courant i,, donnés en vertu de ce qui pré- 
cède, par 
Pin) 


+ RC = =, Lo = — 
Vo Vo 


II. — Discussion de l'équation. 


À. Composantes de l'impulsion. — Soient { — 0 
l’instant de l’arrivée de la particule ionisante dans 
la chambre (sa durée de parcours est supposée négli- 
geable), { — + un instant quelconque où l'impulsion 
est étudiée, { = 7, et { — 7, les instants respectifs 
où disparaissent les derniers ions négatifs et les 


derniers ions positifs, par neutralisation sur les 
électrodes. 

Si pit) et p,(r) sont les parties de la puis- 
sance pr) relatives respectivement à chaque 


catégorie d'ions, l’amplitude u(r) de l'impulsion 


est telle que 


u(r) = ut) u2(T), 


üy et u, étant les solutions de l’équation différentielle | 
lorsque l’on y remplace p(f) successivement par me. 
et Po 

Les résultats généraux obtenus en étudiant & 
seront valables pour &,, et l'étude de na 
totale s’en déduira immédiatement. 


“B. Étude de la fonction Us (29 Solution. 
Avec la convention u, = o à l'instant é— 0 (ui 
libre électrique), PÉte (1 bis) admet la solos 


ess 


Cette fonction présente après r,, instant de dispae 
rition de la dernière charge de la catégorie, une 
décroissance exponentielle du résidu de ten 


sion w(r;), suivant en 
Pas à 


ui(T 5) = (te RE 


T—t 


opte de eo 


ES) Le bis). | 


20 Cas où RC est grand ou très grand par rapport 
à 7. — En posant | 4 


Pape M | 
et en remarquant que e Cest constamment compris 
T 5 F LR 
—-7C 4 
entre r ete © , l'on peut écrire : Re: À 
ne RC OR s 
d’où l’on bre, lorsque RC > 7; : + 
œ ë 
a) = DA — 0), > (01 
avec 5 
S : x2 S De 
PE RO CRC UC ne LS 


L'on en conclut que, lorsque la constante de temps 
du circuit est grande devant la vie +, ge ions, 


l'amplitude de l'impulsion est, à moins de == = près, : 


proportionnelle au travail de ces ions; la Re | 


partie de la fonction (r> 7,) reste alors voisine … 


on) — U, (x), sauf lorsque l’ordre 


de grandeur de + approche celui de RC. 

Lorsque 7, n’est pas très petit devant RC, on peut - 
obtenir le travail w,(r) en intégrant l'équation; il 
vient 


de la constante. 


wi(z) 
CVs 


HA) na : (4) 


RC 


u, (7) étant l'amplitude moyenne entre { — oett —+. 
Dans ce domaine de constantes de temps, nous 
dirons que l'impulsion est du type «travail ». 


SAS PAR EE IN 


A0. Cas où RC est petit devant 7,. — L' impulsion est 
_alors voisine de la fonction PE) comme on s’en 
rend compte immédiatement. sur la figure 2 qui 


#1 
4 


Fig. 2. 


ss interprète graphiquement la relation a bis). L’impul- 
_sion apparaît comme en retard de RC sur la fonc- 


R 
_tion y, P: notons aussi que les maxima et minima 


_de ne sont des points de la courbe Pi 
__ La deuxième partie, 
plus vite que RC est plus petit. 

On serait tenté de réduire la constante de temps 
jusqu’à ce que l'écart entre les deux fonctions 
devienne négligeable ; mais l’amplitüde, en moyenne, 
se trouverait alors réduite dans la même mesure 
que RC: en effet l’ équation (4) ci-dessus, multipliée 


ss RC 
- par - donne 


Ua (T) + RC ts) = 1 Ualr), 


- sauf si l’amplitude instantanée est beaucoup plus 
. grande que l’amplitude moyenne, le deuxième terme 
“du premier membre est petit devant u,(r), de sorte 

_que l’amplitude HS pour RC petit devant 7 
DS Su 1(7); 
devant l’amplitude au 


est de l’ordre de c'est-à-dire petite 


Di(s) =, e 


Comme il est impossible, en raison du bruit de 
fond, de réduire arbitrairement l'amplitude du 


phénomène initial, l'emploi de courtes constantes 


de temps ne donnera pas d'image parfaitement 


_ fidèle de la fonction descriptive du mouvement des 


ions. 

De même que l'impulsion elle-même décrit plus 
-ou moins la puissance, l’aire qu’elle délimite avec 
laxe des temps est liée au travail des ions; en 
intégrant l'équation (1 bis) il vient en effet 


J'utwba+rcuts= 7 wi(r),  (4ter) 


0 


relation qui montre que le travail des ions est propor- 


les fonctions u, ne 


après Ts MO oc 


(Abis) 


_ La 
A 
n B. 
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tionnel à l’aire de l’ impulsion de zéro à +, augmentée 
d’un rectangle de base RC; en particulier si l’on 
Er ou (T0, DC) WT). et"les ‘aires 
totales délimitées ne -dessus de l’axe des temps par : 


x Pi s’équivalent; si À, est l’aire 


. délimitée par l ot nous aurons 


4 
æi(Ti) = FR Ar (5) 


Dans tous les cas où RC est assez petit pour que 


impulsion se rapproche de la fonction p, nous 


dirons que cette impulsion est du type «puissance ». 


4° Cas où RCest du même ordre que *,. — L’impul- 
sion est alors d’un type « mixte », et bénéficie des 


résultats généraux obtenus dans les deux premiers 


cas. Aïnsi en utilisant, soit la relation (4), soit les 
équations (1 bis) ou (4fer), selon que RC sera 
supérieur ou inférieur à 7,, l’on pourra restituer 
commodément les fonctions p,(r) ou w,(r) qui 
décrivent le mouvement des ions. 

Le début de la décroissance exponentielle finale 
sera observable et pourra fournir une mesure appro- 


_chée ou confirmative de RC. 


Quant à l’amplitude de l'impulsion, elle restera 
de l’ ordre de grandeur de 
AE ea ®) pour st AC; 


T CV 


mais elle sera assez sensible à l’allure de p,. En 
particulier si la fonction p, est croissante, u,(r;) 
sera plus élevé que si p, était décroissante. Pour 
fixer les idées supposons RC = 7; ‘la première 
inégalité obtenue au paragraphe 29 ci- ere montre 
que 


Ui(mu) > mu) >= ; Vitm)> ee 


si, par exemple p, était constante, u,(r,) serait égal à 


ET 


Uitn) + 0,6 Ui(ri). 


Enfin, et cela est valable quelle que soit la cons- 
tante de temps, la pente initiale u,(o) est égale 
SP: (0) : 

CV 


50 Influence des paramètres R, C et Vs. Lorsque 
l’on prend ies dérivées partielles par rapport à R 
et C de la solution (2), l’on obtient 


dui(r e se : ue 
je ee r. (x — 1) pi(t) dt, 
( 
dui(s) Ut _ 
ee = C u,(t)eCd 


0 


La première est évidemment toujours positive. 
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Quant à la seconde elle est toujours négative car 


T = T—t 


et 


T T 4n 


l 

outre Ace 7 f u(t) e“C dé. ce 
0 

On a donc toujours avantage à diminuer le plus 

possible la capacité totale; en particulier lorsque RC 

est grand, l’amplitude Ne finale s’identifie 


au quotient bien connu £ (où q représente la charge 


totale des ions) et varie par conséquent en raison 
inverse de la capacité. D'une manière générale on 
s’attache à réduire C le plus possible, et c’est en 
définitive R qui seule commande la constante de 
temps. { 

Quant à V, il n’est cité que pour mémoire car sa 
variation est sans action sur l’amplitude de tension; 
en effet si V, apparaît bien en dénominateur dans 
la solution u,(r), il ne faut pas oublier qu'il est 
implicitement contenu en facteur dans les fonc- 
tions w, et p,. Par contre V, influe sur les durées + 
ou 7,, toutes choses égales d’ailleurs, et, par consé- 
quent, sur le rapport de ces durées avec RC. 


C. Impulsion résultante. — [L'’impulsion par- 
tielle uw, relative à la deuxième catégorie d'ions 
obéit exactement aux mêmes lois que u,, en faisant 
intervenir les fonctions p,(r) et w,(t), ainsi que le 
paramètre x,, instant de disparition du dernier ion 
de la catégorie. 

L'impulsion résultante u, somme de u, et w, 
obéit également aux lois déjà énoncées pour w;; 
il suffit de remplacer partout u, par u et p, par p. 

On peut distinguer dans le développement de 
l'impulsion totale trois phases principales; sir; < 7%, 
on rencontrera : 


Une phase initiale où o <7r <7 


DD L 


u(r)= —— ‘ eRC[ pi(t) + pa(t)] dé 
CA 
Une phase intermédiaire de 7, à 7, : 


(n 


CRAN CRô 
ref * pa(t) dt; 


Une phase finale de 7, à l'infini 


ufr) = tn dt de 


LR) = en 


Par des moyens appropriés (collection d'électrons) 
l’on peut obtenir pour les ions négatifs des vitesses 
de déplacement 10% fois supérieures à celles des 
ions positifs, ce qui permet de séparer très aisément 
dans le temps les influences respectives des ions 
de chaque signe sur l’impulsion. 
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D. Impulsion complémentaire. — En suppo- 
sant les paires d'ions distribuées exactement de la s 
même manière à l'instant { — o, et si l’on inverse 
la tension V,, on obtient une impulsion négative 


égale à 
T 


te 
——— 
CVs à 
dæ; d#o 
ais ; 
dre dé 


t 
ui(r) =— eXC p;(t) dt, 


avec 


PIS PUREPRES 
les durées 7; et x, de vie des ions seront alors 
différentes de 7, et x, comme pour la puissance ef 
le travail. 

Cette impulsion. négative , cependant liée Le 
l’autre par les relations 


Wa(tu) =#w» (7), 


wi (ti) ee Wjo(Tio), 


Wat 


qui résultent du fait que les parcours et chutes de. 
potentiel des ions sont échangés entre les deux 
catégories quand on inverse la tension, en négligeant 
naturellement l’écartement initial des deux ions de 
chaque paire, effet qui d’ailleurs se compense statis= 
tiquement. : 


IT. — Amplification de l'impulsion. 


Par un raisonnement tout à fait analogue à celui 
qui nous a permis d'établir l'équation différen- 
tielle (1), on peut étudier sans difficulté la transfor-. 
mation subie par l'impulsion lorsqu'elle traverse 
une liaison résistance-capacité de l’amplificateur. 
Si u; est la tension transitoire reçue par la grille 
de la lampe n°0 i, — k; u; la variation correspondante 
sur la plaque de la lampe, l'impulsion u;,, que reçoit. 
la grille de la lampe suivante est donnée par 


! 
= — K;u;, 


R;41 Cia = 0,1 étant la constante de temps du COU= 
plage et k; le coefficient dynamique d’amplifica- 
tion de la lampe n° i au point de fonctionnement 
considéré. | 
D'après la discussion déjà faite pour l’équa- 
tion (1 bis) il apparaît que si 0;,1 est très grande 
devant la durée du phénomène transitoire, l’on 
aura une amplification pure et sans distorsion 


= — kiu:. 


ui 


Si, au contraire, 0;,, est du même ordre que la durée 
de l'impulsion u; ou d’un ordre très petit, w,1 
suivra plus ou moins la dérivée u/, avec un affaiblis- 
sement d’ amplitude du même ordre que celui du 
rapport de 0;,1, à la durée de l’impuision; en outre 
une distorsion de phase comportant un retard égal 
à 0;,, interviendra. 


Pad 5e è : 2 
a" f SA È 
va ” - 4 + y 


_ C'est pourquoi généraiement l’on recherche la 
fidélité en maintenant les constantes de temps de 
tous les étages, y compris celle de la chambre d’ioni- 
sation, assez grandes par rapport à la durée de 
l'impulsion, sauf une dans certains cas. 
Supposons que l’on ait opéré ainsi et que 0 soit 
cette courte constante de temps. Dans ces condi- 
tions la tension transitoire de sortie de l’amplifi- 
-cateur est en valeur absolue 


T 
Le u 


(ere E 
Un) = KT [ Pt) el dé, 
20 


K étant le coefficient total d'amplification, compte 
tenu du feed-back s’il y a lieu. 
u,(T) est en somme la solution de l’équation 


_— + U, = KP) ({ter) 


Ch 


qui se discute exactement comme (1 bis), u, rempla- 


-çant u et Ü remplaçant RC. 


Remarque. — Dans une liaison résistance-capacité 
entre deux étages d’un amplificateur, la résistance 
et la capacité interviennent inséparablement par leur 
produit, donnant ainsi un degré de liberté dans la 
fixation de la constante de temps. C’est une diffé- 
rence importante avec le circuit initial (chambre 
_d’ionisation) où la capacité C influe directement 
sur l’amplitude. 


IV. — Exploitation de la théorie. 


- Cas des géométries à champ droit. 
- À. REMARQUE PRÉLIMINAIRE. — Les lois énoncées 
dans ce qui précède sont valables pour une géométrie 
quelconque, quelle que soit l’histoire des ions. Ce 
résèltat n’a pu être atteint précisément qu’en 
incluant implicitement ces importantes conditions 
dans les fonctions p et w. Si l’on veut continuer 
l'étude des impulsions en explicitant certains carac- 
_tères de ces dernières fonctions. il devient indispen- 
_sable de faire intervenir la géométrie de l’appareil, 
d’une part, et, d’autre part, l’ensemble des lois 
régissant le comportement des ions dans les gaz 
(mobilité, recombinaison, attachement, diffusion, 
jionisation par chocs, etc.). Cette dernière partie 
nécessiterait à elle seule une étude particulière que 
nous ne pouvons entreprendre ici. Nous suppo- 
serons que l’on a. réalisé des conditions expéri- 
mentales telles que le nombre des ions libérés par 
une particule ionisante ne soit pas affecté de manière 
appréciable par l’attachement et la recombinaison 
et qu'il n’y ait pas d’ionisation par choc. 

Nous allons montrer ci-après quels renseignements 
sur la particule ionisante il est possible de tirer 
de l’étude d’une impulsion. 
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B. DÉTERMINATION DU NOMBRE DE PAIRES. 
On sait que le nombre de paires est lié étroitement 
à la perte d'énergie de la particule pendant son 
passage dans le volume gazeux soumis au champ 
électrique de la chambre, d’où la possibilité d'utiliser 
un tel appareil comme spectromètre. 


19 Loi générale. — Le travail total des ions créés 
par une particule dans une chambre est propor- 
tionnel au nombre de paires. Les méthodes de 
restitution de la fonction w(r) envisagées précé- 


. demment permettent donc de déterminer dans tous 


Jes cas le nombre de paires à un facteur près, si 
l’on possède un enregistrement convenable de 
l'impulsion issue de la chambre. 


- 20 Arlifices permis par la collection des électrons. 
— Lorsque l’on utilise la méthode de collection 
des électrons et que l’on élimine par un choix judi- 
cieux de la constante 0 l’influence des ions positifs, 
il est encore possible de déterminer dans certains 
cas le nombre de paires. 

Ainsi Sherr et Peterson (?) rappellent ou suggèrent 
quelques artifices applicables le plus souvent à la 
seule géométrie plane (électrodes planes parallèles). 


Emploi d’une grille (3) (applicable à d’autres 
géométries); 

Emploi d’une électrode fine (inapplicable à la 
géométrie plane); 


Méthode de la différentiation de l'impulsion; 
Attachement retardé des électrons. 


Il faudrait y ajouter l’utilisation d’un champ 
électrique alternatif. Ces auteurs ont employé avec 
succès la troisième méthode qui consiste, en accord 
avec nos propres conclusions, à obtenir une impul- 
sion du type puissance, dont le palier initial mesure 
par son ordonnée le nombre de paires; ‘ils tentent 
également d’étudier lorientation des trajectoires 
par rapport au champ et l’ionisation spécifique 
des particules. Soulignons que cetté méthode est 
limitée strictement au cas de la géométrie plane, 
d’ailleurs la plus simple et la meilleure dans la 
plupart des cas. 


C. ORIENTATION DE LA TRAJECTOIRE PAR RAPPORT 
AU CHAMP LORSQUE LES LIGNES DE FORCE SONT 
DROITES. — 19 Méthode des discontinuilés. — Nous 
allons envisager le cas où les électrodes sont, soit 
des plans parallèles, soit des cylindres concentriques, 
soit des sphères concentriques. Les lignes de forces 
sont alors des droites, et les surfaces équipoten- 
tielles des plans, cylindres ou sphères parallèles 
aux électrodes (en mettant de côté bien entendu 
les régions extrêmes où les bords des électrodes 
introduisent une distorsion). 


(2) Voir références bibliographiques. 
(#) ALrEN, voir références bibliographiques. 


j- 
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Dans ce cas, toutes choses égales d’ ailleurs, les 
ions de chaque catégorie prennent un mouvement 
moyen bien déterminé qui peut toujours se traduire 
par une relation théorique ou expérimentale entre 
la position sur la ligne de force et le temps; dans 
des conditions invariables (gaz, pression, tension) 
un ion de telle catégorie parcourra de la même façon 
en fonction du temps la ligne de force sur laquelle 
il s’est trouvé créé, quelle que soit cette ligne de 
force, aux fluctuations d’agitation thermique près. 

D'autre part la fonction p(r) relative aux deux 
catégories d’ions, changera en général quatre fois 
après l'instant initial : 

Aux instants f, et { où les ions de chaque catégorie 
commencent à disparaître par neutralisation sur 
leurs électrodes respectives; 


Aux instants z, et Tr où les derniers ions de 
chaque catégorie disparaissent et après lesquels les 
fonctions p, et p, restent nulles. 


L'on conçoit que si l’on peut repérer sur l’impul- 
sion ces instants, on puisse connaître sur quelles 
surfaces équipotentielles se trouvaient initialement 
les deux extrémités de la ligne des paires; ce rensei- 
gnement, joint à la longueur de la trace déduite 
du nombre de paires et de la relation énergie- 
parcours ou plus exactement ionisation-parcours 
relative à la particule dans l’atmosphère gazeuse 
considérée, nous donnera, par exemple, par cons- 
truction géométrique : 


L’angle de la trajectoire avec le champ si la 
géométrie est plane ou sphérique; 


Une relation géométrique entre deux paramètres 
d'orientation de la trajectoire par rapport au champ 
et à l’axe des cylindres, dans le cas de la géométrie 
cylindrique. 


* 


Mais les changements d’allure de la puissance ne 
donnent sur l’impulsion des points anguleux (les 
plus faciles à déceler) que dans le cas expérimen- 
talement difficile à atteindre où cette impulsion suit 
suffisamment bien la puissance; en pratique ces 
points anguleux seront d’ailleurs toujours arrondis 
par la diffusion statistique des ions autour de leur 
mouvement moyen théorique. En fait les rensei- 
gnements obtenus par cette méthode seront géné- 
ralement très imprécis. 

Le problème se simplifie considérablement lorsque 
l’on connaît déjà la position en potentiel de l’origine 
de la trace, en particulier lorsque la particule est 
issue de l’une des électrodes. Nous allons indiquer 
ci-après une méthode valable dans ce cas, et qui 
nous semble susceptible d’une PRÉCEIQR bien supé- 
rieure. 


20 Méthode du centre des charges. — Supposons 
que l’on ait pu déterminer, en plus du travaii total 
des ions w(t;), le travail partiel d’une des caté- 


_gories d'ions, par exemple W; G) . ions négatif 

Si r est le paramètre de longueur qui perme 
repérer la position d’un ion sur une ligne de | 
(distance à l’une des électrodes planes, ow à l'axe 
des cylindres, ou au centre des sphères), on sait. 


que le potentiel V d’une surface équipoten- | 
tielle r — const., peut se mettre sous la ce | 
suivante : 3 


Le = A (T5 — r), pour la sédréteie plane ; 
— a (Log r, — Logr) pour la géométrie sine. 
Hi. 
I 


V=a Î — =) pour la géométrie sphérique; | 


étant la valeur que prend r pour l’électrode a 
potentiel zéro, et a un coefficient qui dépend de. la. 
tension V, et de la dimension de l’autre électrodes 
Une charge quasi ponctuelle Ag; se trouvant. 
initialement sur la surface équipotentielle F SP 
effectuera un travail total 


a(ro— pi) Ag TER 
ou ne 
a(Logro— Logp;)Agi  - 
ou 3 a Wu > 
«(= = =) Ag;, si elle est positive . ie 
Pi To / EN 
7 
ou la différence de ce travail avec Vo Ag: si el E 
est négative. Es 
Le travail w,(r) de toutes les charges positive S. 
sera donc, selon la géométrie, €. 
Æ 
alto) —:a ro—©;)Ag; = a —_à : Ag;  : 
(2) > o— pi)Aqi= arog —< 2° 1° 2 
na > ex 
ou { 
wa (52) = à Ÿ (Logro— Logps) Agi % 
| #. 
= a Log ro q — D Loge; Ag: E : 
î L 


© 
[er] 
Fr 


e(m)=a) (EE )ansa taf; F- 
L 
i i 


.q étant la charge totale ZX Ag; proportionnelle au 


nombre de paires. On aurait par différence le travail 
des charges négatives w,(7,). 
Soit p la grandeur définie par 


De: Agi 


l » , ® 
ee Fa … pour la géométrie plane; 


2e 


L 


ce pour la géométrie cylindrique; 


_ Dee I 


(”) Voir l'exemple traité ci-après. 


pour la géométrie sphérique. 


trajectoire un point que nous appellerons « centre 
des charges » et tel que le rapport des travaux 
partiels des ions soit le même que s'ils étaient tous 
créés initialement en ce point. Ce point se confond 
avec le centre de gravité dans le cas de la géométrie 
plane. 

La variable 9 est en général fonction : 


De la particule (relation ionisation- -parcours dans 
1 atmosphère considérée) ; 


_ De l'énergie de la particule; 
: WRe : RO de l’origine de É trajet: 


Elle est indépendante de la tension V,. 
… Laissons de côté le cas du cylindre trop complexe 
et à peu près inutilisable lorsque le cylindre intérieur 
ne se réduit pas à un fil. 

À. Il est facile, en s’aidant de ce que l’on sait par 
ailleurs sur l’ionisation des particules, et au besoin 


par des expériences spéciales préliminaires, de déter- 


dentiques d'énergie constante issues d’un point 
xe pris sur une électrode. 

_ Si l’on étudie une impulsion produite par une 
articule de type connu ou identifiable, si l’on est 
certain que l’origine de sa trace se trouve sur une 
électrode (ce que pourra confirmer la position de la 
discontinuité correspondante), et si l'on a pu déter- 
miner w et w,, par conséquent =? ——"®-, on 

wi Œ — Wo 

“connaîtra aussi le nombre total de paires (donc 
l'énergie), la- surface lieu des centres des traces 
des particules de cette énergie, et par l'intersection 
de cette surface avec 1a surface équipotentielle 


‘dont le paramètre r correspond au RATE . 
1 


mesuré, l'angle © cherché. 

Cette méthode est inairecte, plus compliquée et 
plus restrictive (nécessité de situer au préalable 
Vorigine de la trace) que celle des discontinuités. 
Elle est toutefois très aisément exploitable dans le 
cas de la géométrie plane, et susceptible de donner 
Le bons résultats. 


ke 
- D. EXEMPLE D'UNE CHAMBRE A  ÉLECTRODES 
PLANES PARALLÈLES. — On se propose par exemple 
d'étudier les impulsions données par les particules 
alphas émises par un dépôt mince sur la région 
centrale d’une électrode. La chambre est disposée 
et alimentée en tension comme J’indique le schéma 
de la figure 5. 

_ La pression est telle que le parcours des particules 
Bu inférieur à la distance des électrodes. Le gaz 
utilisé est tel que les électrons restent libres. 


. 


| IMPULSIONS DE TENSION ISSUES DES CHAMBRES D’ LE 


Cette lee 5 du paramètre r détermine sur la - 


miner la surface lieu des centres pour des particules 
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La tension est assez grande pour assurer la satu- 
ration ; nous négligerons la recombinaison résiduelle. 
L'électrode collectrice est reliée à un amplifi- 
cateur dont la constante de temps 0 est du même 
ordre que le temps de collection x, des ions positifs. 


Fig. à. 


Les impulsions de sortie sont appliquées aux 
plaques verticales d’un oscillographe cathodique; le 
balayage horizontal est déclenché par un dispo- 
sitif approprié par l’arrivée de l’impulsion elle-même, 
et sa période est choisie de l’ordre de 1 à Ta; On 


photographie les impulsions. 


En faisant abstraction : de 
_ De la diffusion des ions; 
- De leur recombinaison ; 

Du bruit de fond ; 


_ le calcul permet de prévoir des courbes 1, 2, 3 (fig. 4) 


Fig. 4. 


relatives aux traces 1, 2, 3 de la figure 3, ramenées 
au même instant initial. Re 


19 Constatations immédiates. de Les courbes 2 et 3 


: présentent un décrochement initial dû au mouvement 


des électrons; ce mouvement est très rapide relati- 
vement à la vitesse du balayage.et apparaît comme 
instantané. Du fait que 0 se trouve très grand 
devant la durée 7, de vie des électrons, la hauteur 
du décrochement est proportionnelle au travail 
des électrons w;(Tt:). 

La courbe 1 comprend seulement deux branches 


! 
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exponentielles séparées par l'instant *. tandis que 
les deux autres présentent une branche intermé- 
diaire allant de £ à 7. 

Les discontinuités en {, et même en 7, (sauf pour 1) 
sont presque indécelables, bien que l’on n'ait pas 
tenu compte de la diffusion; les discontinuités en°To 
sont dans le cas présent des points d'inflexion pour 
les courbes 2 et 3. 

Par contre, les aires apparaissent tout à fait 
mesurables; le bruit de fond, en épaississant le trait, 
ne pourrait, sauf s’il était exagéré, entacher cette 
mesure d’une erreur très importante, alors qu'il 
interdirait complètement de repérer les discon- 
tinuités. 


20 Mesure des énergies. — On mesure les aires 
totales depuis l'instant zéro jusqu'à l'infini; comme 
on ne peut observer, en raison de la période du 
balayage, la totalité de la courbe, on remplacera 
la partie manquante par un rectangle de largeur 
égale à 0. 

Notons que, si l’on est certain que les particules 1 
et 2 sont d’égales énergies, il est possible de déter- 
miner 0 graphiquement en cherchant la largeur du 
rectangle à ajouter à la partie de droite de l’are 
comprise entre les deux courbes, de façon que l’aire 
ainsi obtenue égale la partie de gauche. 


V 2? 


Si À est l’aire totale de l’une des courbes (fig. 5), 
le travail total des ions est 
| CM 
0 


A. 


Mi Wii + Po = 


On en déduit le nombre de paires N total 
Ne 


| 
€ Vo di) è 


où € représente la charge élémentaire; du nombre 
de paires il est facile de passer à l’énergie. 
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L Li LL ner "és "2, 7 ns! 2. ho 1 
v79 ñ LS LT DEA Ne ATOS + | 
s Az \ ; y ILE PONS | 
\ « : o ; rh LISTER + 
’ ‘ ; NH À x 2 
+ Fe CIE 


30 Détermination de l'orientation de la trace (fig. 5). 
— [Le décrochement initial w, représente, compte 
(Ti) 


tenu de l’échelle, la quantité . On peut donc 


CV 
déterminer le quotient é 
à wi(T1) LED 120% , 
W1(T1) + W2(T2) SCA 14 


Par ailleurs, la courbe de Bragg relative au gaz 
de la chambre, et modifiée pour tenir compte de la 
pression, nous permet de situer le centre de gravité 
des charges par rapport aux extrémités de la trace, 
cette position relative variant de façon déterminée 


_ avec l’énergie de la particule. Ayant déjà déterminé 


l’énergie des particules étudiées, nous pourrons alors | 

connaître la distance / de ce centre à l’origine de la. 

trace, et, par conséquent, l’angle « qu’elle fait avec 

le champ sr 
d wo 


Vibes PAT 


d'étant l’écartement des électrodes. de ” 
On. remarquera que cette opération n’a fait inter= 
venir que les travaux totaux des ions de chaque 
catégorie, et par conséquent, qu’elle doit rester 
insensible à une diffusion des ions, même impor- 
tante, ce qui n’est pas le cas pour les discontinuités. 


4° Autres renseignements. — On pourrait se. 
livrer sur les courbes présentées, à la restitution. 
complète de la fonction p,;, puissance des ions 
positifs, en ajoutant à l'amplitude en chaque point 


_ un supplément 0 obtenu graphiquement (fig. 5). 


La courbe ainsi obtenue représenterait, comme l’ont 
indiqué Sherr et Petterson, l’ionisation en fonction 
du parcours, dans sa partie décroissante. 2 
Nous ne pensons pas que cette opération, dans le 
cas présent, puisse être assez précise pour fournir 
autre chose qu’une allure ou un ordre de grandeur. 
Signalons, en outre, en ce qui concerne la recombi- 
naison, que l’effet de colonne devrait croître quand 
l’angle © diminue, et que les méthodes suggérées 
ci-dessus pourraient être appliquées à l’étude de ce 
phénomène. : 


V. — Conclusion. 


« 


À la lumière de la théorie sommaire ainsi élaborée 
il n’est pas interdit d'envisager que les chambres 
d’ionisation, moyennant peut-être certains -perfec- 
tionnements, puissent être utilisées comme de bons 
spectromètres, et même comme des appareils 
descriptifs de particules. Dans ce dernier ordre 
d'idées, les renseignements obtenus resteront sans 
doute limités et pour certains, imprécis; en revanche, 
dans tous les cas où J’on peut s'en contenter, la 


RÉSISTANCE EN 


: 


Hombre d’ionisation présenterait d’incontestables 
avantages de rendement et se prêterait à des 
recherches statistiques sur des phénomènes peu 
fréquents, à l'encontre d'appareils beaucoup plus 
complets et plus précis comme la chambre de Wilson. 
Cette technique peut s’accommoder de circuits 
électroniques relativement simples donc peu encom- 
brants et semble susceptible de rendre des services 
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dans le domaine du rayonnement cosmique en 


particulier. 


Je tiens à remercier ici M. G. Ambrosino, dont 
les suggestions et les critiques avisées m'ont conduit 
à préciser certains points importants de ce travail, 
ainsi que M. Maurice de Broglie pour ses précieux 
encouragements. 


Manuscrit reçu le 28 avril 1949. 
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ÉTUDE DE LA RÉSISTANCE EN HAUTE FRÉQUENCE D'UN ENROULEMENT A FIL DIVISÉ 


Par A. COLOMBANI. 
Faculté des Sciences de Dijon. 


Sommaire. — Après l’étude magistrale de Sommerfeld concernant la résistance en haute fréquence 


. des formules approchées s 


d’une couche cylindrique et uniforme de matière conductrice, Butterworth puis Austin, 
s'appliquant à des enroulements à fil plein ou divisé en brins. 


établirent 


En partant directement des équations de Maxwell et en supposant le rayon du « brin » inférieur à 
l'épaisseur de pénétration du courant, j'ai pu établir une formule très simple concordant parfaitement 


avec les calculs de Sommerfeld et avec l'expérience. 


: En l’appliquant à différents problèmes il m’a été possible de mettre en évidence pour chacun d'eux 
*. l'existence d’une résistance optima et d’une forme particulière de bobine. 


Les résultats théoriques obtenus concordent avec les données pratiques fournies par la technique 


moderne des enroulements. Il m'a paru intéressant de les signaler aux physiciens qui ont à déterminer 
rapidement et avec précision des bobinages répondant à des conditions données. 


_. Considérons un conducteur sensiblement rectiligne 
et éloigné des conducteurs voisins parcouru par un 


courant de haute fréquence. On sait que la distri- 


bution du courant dans sa section droite n’est pas 
uniforme. Plus la fréquence est élevée, plus la 
densité de courant augmente en allant de l’axe vers 
a surface. 

- Cependant la répartition du courant autour de 
l’axe du conducteur est régulière et l’accroissement 
de résistance qui en résulte est aisément calculable. 
En particulier si le rayon r du conducteur est petit 

I 
Vis-à-vis de « l’épaisseur de peau » € — ere 
résistance calculée est sensiblement la même qu'en 


courant continu et la répartition du courant est 
uniforme dans la section du conducteur. 
Envisageons à présent le cas de fils rapprochés 
parcourus par un courant de haute fréquence. Le 
champ magnétique produit à l’intérieur d’un fil 
par les courants de haute fréquence voisins cesse 


d’être négligeable par rapport à celui que le fil lui- 


même produit dans son intérieur. La distribution 
du courant dans le fil est donc très différente et la 
résistance apparente plus grande que ne le ferait pré- 
voir l'application des formules de l'effet Kelvin 
relative à des conducteurs très éloignés. 

C'est ce qui se produit dans tous les bobinages 
et en particulier dans ceux dont le fil est « divisé » 
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afin d'utiliser au mieux le métal parcouru par le 


courant de haute fréquence. Cette « division » est 


réalisée au moyen de brins très fins, isolés et 
« toronnés » ou « cordés ». Leur rayon individuel 1k 
satisfait à la relation r<<: et leur enroulement 
est tel que chacun- d'eux occupe successivement 
dans la section du fil toutes les positions possibles. 
De cette facon ils sont tous placés dans les mêmes 
conditions vis-à-vis du champ magnétique global 
extérieur, Et l’on peut très sûrement admettre que 
le champ magnétique dans chaque brin est princi- 
palement dû au voisinage des autres brins du même 
fil et à la proximité des spires juxtaposées du bobi- 
nage qui ajoutent leurs actions à l’intérieur du brin 
considéré. | 

Autrement dit la condition d’uniformité du champ 
dans lequel le brin est placé se trouve réalisée avec 
une grande approximation lorsque le champ magné- 
tique dû au brin lui-même dans son intérieur est 


petit par rapport au champ H dû à tous les autres 
brins du bobinage. 

Supposons donc un champ magnétique alternatif 
d'amplitude constante H perpendiculaire à l’axe 
du brin (fig. 1). (Nous verrons plus loin l'influence 


Fig. 


de l'obliquité du brin sur le champ). Nous pou- 
vons assimiler à un conducteur rectiligne une 
portion de brin suffisamment longue par rapport à 
son diamètre car celui-ci est très petit vis-à-vis 
du rayon moyen d’enroulement de la bobine. Prenons 
dans la section droite du brin deux axes rectangu- 
laires O:, O; l’axe O, étant dirigé dans le sens du 


champ H. L’axe du fil portera O, parallèle au champ 


e ve . 57 ‘ 
électrique À suivant lequel s'écoule le courant 
électrique. 


+. < 


NO er te). | Re 
La relation roth — donne immédiatement 


oh 
RS o 11. 


Soient 1, la densité de courant et y la conducti- 


Tes 


bilité. On a. LE 


D'où avec H uniforme 
ta = td, + ©Y Hz. 


Par intégration sur la surface circulaire qui cons- 
titue la section droite, on trouve que l'intensité 
totale dans le brin de rayon r est i = tri, D'où 


: Ü ? He 
la =: et par conséquent Re 


7? 
: # > “ à 4 
= +u Hz: RE  - 
zr EC. 4 


Ce résultat est d’ailleurs évident car la valeur 
moyenne de z est nulle et t,, est la densité moyenne 
de courant dans la section droite. 

On en déduit l’effet Joule par unité de longueur 


: Î 15 ds, ds étant l’élément de surface 
dS —— 2R? sin? da avec. 3 —7rCos&r: 
Soit, tous calculs faits 


12 ro?yrt H? 
+ ——— + 


Le premier terme En est l’effet Joule normal 


produit par l'écoulement du courant à à travers 
la surface Tr? sur une longueur unité. | 
Le rapport de la résistance en haute fréquence 


à la résistance ohmique a donc pour valeur 2" 
ARE 7? Ææ DS 
< 700) Ve x 
TES M TE .@ 


H? est la valeur moyenne du carré du champ magné- 
tique dans le volume occupé par tous les brins de 
l’enroulement. L’intensité à est supposée la même 
pour tous les brins et constante tout le long de l’en- 
roulement. Cela nécessite donc l’utilisation de la 
bobine sur des fréquences très différentes de sa fré- 
quence propre afin d'éviter la production d’oscil- 
lations localisées et par conséquent de variations 
d'intensité le long des fils. Enfin comme nous le 
verrons plus loin il ÿ a intérêt à ce que la capacité 


répartie soit rejetée aux extrémités de l’enrou- 
lement (1. 


ES. 


Afin de calculer le rapport # nous allons main= | 


12 


tenant avec Paul Langevin faire intervenir la forme 
de la bobine à l’aide de trois coefficients : &, Gare pd. 

Supposons que l’enroulement comporte n spires | 
formées chacune de N brins. Soit D son diamètre 
moyen. La section droite de la “eo de la bobine 


(1) Quand les fréquences sont voisines et la capacité rejetée 
aux extrémités on prend le point milieu de l’enroulement 
pour mesurer l'intensité qui servira à définir la résistance 
pour la fréquence fondamentale. On pose W = R1, W étant 
la puissance totale. 


Le volume correspondant à cette gorge est 


V=zxSD= Tropi. 


: RTE TT 
ne one 


'A'A'A'A'A'A'A'A'A' 


nt 
qe QU 


il 
AN 


- DRE TN 


D'autre part la résistance et la self de la bobine 
en courant continu ont respectivement pour valeur 


D 


et L=8nr2D 


L'énergie magnétique localisée dans la gorge est 


LH 
3 D: 


I rs 
MD Une 
a 


_ L'énergie magnétique te est 


- avec Fe Nu. 


« Appelons n leur rapport plus petit que l'unité. 
Ona a 


IP Ds = My e (3). 


I & 


es. = : 

Les coefficients x, 6, n ne dépendent que de la forme 
de la bobine et non de ses dimensions. 

- On en tire donc 


(3) 


- (2?) Avec une gorge circulaire de rayon R on aurait 


TR? mr? D 
LL —) ES 


D? 2R. 


(5) Nous poserons = _. L'expérience justifie cette hypo- 


thèse (voir plus loin). 
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= Par conséquent (2) devient 


n D n ZI? 
Rp—= Nr Li x? r0 «2? Le TUE | 
RDS F n?/V? 
= 2 7602 y2 IE 
ne [ir PONTS Ê 
ou encore 
G/ x 2 
hr Roi F (5 rwrnN) | (4) 


ARE I ue Re S 
Rr= Er. A (5 ré vo). AI (5) 


Remarque. — Considérons la valeur 


Rr= Bof i+ A (gro nN) |. 


Avec les hypothèses précédentes, on a 


LE 
c Re l2 
2 ‘ l ; 
D'autre part prenons Nn — _ ce qui correspond 


à une gorge remplie complètement par le métal de 
l’enroulement, et posons e — 2 gr, q étant le nombre 
de couches supposées jointives. 


On en tire: 


ri 


Rp= Roli+ 47r2y2w?2g2r1] où = fi+at |, 


en désignant par € l'épaisseur de pénétration 
1 


E — (270ÿ) ? 


On retrouve ainsi facilement un résultat établi 
autrefois par Somrerfeld [2] pour un cas A THRIMES 
bien vérifié par 1’expérience. 


Discussion sur la validité des résultats. — 


0 2 : , » C2 2 
10 Considérons d’abord le champ magnétique H. 
Nous l’avons supposé uniforme. En fait, ce n’est 


= 
pas tout à fait le cas. Le défaut d’uniformité sur H 


provient surtout de la partie de ce champ due à 
l’action du courant passant dans le brin lui-même 
et les brins voisins, principalement de ceux qui 
sont toronnés avec lui dans le même fil, L'influence 


. maxima de ces fluctuations sur le résultat c’est-à- 
dire l'erreur relative sur le second terme de (2) 


peut se calculer en remarquant que les variations 
d'intensité du champ dues aux courants dans des 


ere =. , 0 
brins voisins, sont au plus de l’ordre de se 


Comme 


AS 
Ve fn L 


on a donc pour l'erreur relative sur VA? 


2 L )3 
pare VE (U). (6) 
2 Vir IN7 1n 2 
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Cette quantité qui est la condition d'uniformité (AL) 
doit être petite (. 

D'autre part le champ rt serait normal au plan 
des spires si celles-ci étaient centrées sur l'axe de la 
bobine. Mais une composante longitudinale du champ 
magnétique s’introduit par suite de la forme héli- 
coïdale des enroulements. 

Deux effets en résultent : 


A. Le que n’est pas normal au plan ner 


contenant il. 
I en rue par conssauens une composante 


longitudinale , du champ H qui donne naissance 
par induction dans la masse du brin à des courants 
circulaires autour de son axe. 


Les équations de Maxwell s’écrivent pour ce cas 


(fig. 3) : 


Jah (ja= densité de courant), 
1 (rh 7 
rOtAy—= - ( 1): — 11, d’où === y0 
FENOr 2 
et 
: Or 
Ja=— H, (Hy constant), 


r désignant la distance d’un point P de la section 
droite à l’axe du brin. 
Supposons la longueur d’onde de l’oscillation 


créant Îl grande devant les dimensions du conduc- 
teur et négligeons les courants de déplacement, 
Un calcul classique donne pour la puissance dépensée 
par effet Joule [4 : 


pot 


l D 
8 TYW?r* 


2 


(7) 
« * Le 

soit avec H;— H sin0(6angle de H avec la nor- 

male au brin) : 

Ty? rl 


Wi= © 


HP sin? 0, 


r désignant ici le rayon du brin. 


(*) Le calcul de la résistance optima nous en donnera une 
signification très simple. 
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L’angle 0 fait intervenir le pas de l Ro 
Lorsque les spires sont d’une inclinaison telle que sin2ÿ 
n’est pas négligeable, il y a évidemment lieu d'en 


tenir compte dans l'expression générale de Rz 


L'expression (2) devient 


Rp He ro (1 = sin® 1) . | 


Pin Bof + nÊ (2 Tps a (+ : sin: |, 
œ \ D 5 


Y£ I : DER ( ER. )| = 
Es Fa 2 20 Ye 
Br É + 7? NE rSynN|r+ = sin 


ou 


Rr= 


Le 


En introduisant le pas de l’enroulement p —7TDtg0 


on obtient 


LR : 
ie nl. nN RAA 


Si 0 est petit le terme ; Sin? 20 est négligeable CO). 


B. Les N brins du fil enroulés suivant une hélice 
dont la tangente fait un angle © avec l’axe du fil 


produisent à l’intérieur de celui-ci, suivant son axe 


un champ H,— En o, a étant le rayon du fil 


(fig. 4). Ce champ longitudinal donne lieu à un 


dégagement de chaleur calculable par (7) en rem- 


plaçant A, par H, 


ryw?rt {N2r? 


Ws= _ 


sin?0. (9) 


Son rapport avec le second terme de (2) est 


2/V27 . à 
— sin?9, 


a? IT? 


(5) Nous négligerons 0 pour les applications. 


RS | 


LL LES | 
PAL sin? €. (10) 
Comme d> Nr?,W,est supérieur à —— SN sin?o 


&. nn 2 4nL 

_ est-à-dire au carré dé la condition (U) multiplié 
- par 2 Nsin?o. On peut donc HUE que W, est 
a. certainement négligeable. 


20 Nous avons admis que l'intensité t était cons- 


‘4 oo che. la bobine fermée sur une capacité T 
dans un circuit résonant. La force électromotrice 
_ induite dans l’enroulement produit des ondes qui se 
- propagent dans les fils et se réfléchissent partiellement 
- à la discontinuité du condensateur. Il en résulte 
un système d’ondes stationnaires presque parfait 
en raison du peu d’énergie dépensée, 
A la résonance tant que Fest assez grand la lon- 
 gueur d’onde À, le long de l’axe de la bobine est 
grande devant Ja demi-longueur du fil ® et le calcul 


long du fil jusqu’à ses extrémités. Si F diminue, 
la fréquence de résonance augmente, À, 
et tend vers ® et le courant à subit des variations 


dérablement aux extrémités de l’enroulement, et 
. présente un ventre en son milieu. Si nous désignons 


- que possible et supposée rejetée à ses extrémités, 
- la nécessité de maintenir l'intensité constante le 
long de l’enroulement conduit donc à l'utiliser sur 
; des fréquences très différentes de sa fréquence propre 


De 
27 LE 


h= 


En appelant T la capacité en circuit, la condition 
de bonne résonance sera donc T >cC. 


Remarques sur la capacité répartie C. — 
“La capacité répartie doit être aussi faible que pos- 
ble. Une capacité trop élevée limite en effet infé- 
“rieurement la bande de longueur d'onde recevable 
avec la bobine et un condensateur donné T. Elle se 
traduit d’autre part par des pertes d'énergie, c’est- 
“à-dire par une augmentation de la résistance Ry, 
donc aussi de l’amortissement. 
Pour des raisons d’encombrement, les enroule- 
… ments sont réalisés en plusieurs couches bobinées 
“de différentes manières. Souvent les bobinages sont 
- en sens inverse d’une couche à l’autre : la tension 
tre spires correspondantes est done constante 
ce qui diminue les pertes. D'ailleurs un intervalle 


_ tante le long de |? enroulement. Ceci conduit è à consi- . 


montre que l'intensité est sensiblement constante 
diminue 


sinusoïdales le long du conducteur. Il décroît consi-, 


par C la capacité propre de la bobine aussi faible 
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. d’air d’autant plus faible que la fréquence est plus 


basse sépare spires et couches (6). 

En assimilant les couches d’un tel enroulement 
aux armatures de condensateurs cylindriques 
coaxiaux montés en cascade on obtient comme 
ordre de grandeur de la capacité propre 


diamètre moyen de l’enroulement ; 
: longueur axiale; 

: épaisseur de l’enroulement; 

: rapport de l’épaisseur totale de l'isolant (air 
compris) à l’épaisseur totale de l’enroulement. 


ARSAEUCS 


En particulier si l'épaisseur de l’isolant séparant 


I 
reset 
2 


deux couches égale le diamètre du fil 
pour une gorge carrée : | —e. Donc C est de l’ordre 
de grandeur du diamètre moyen D exprimé en 
centimètres. En U. E. M. C est de l’ordre de gran- 
deur de 
=D. (11) 
Il y a donc intérêt à employer des bobines de 
faible diamètre moyen si l’on désire diminuer la 
capacité répartie. De plus l’étalement extérieur 
du champ de la bobine est moins grand avec un 
petit diamètre. On diminue également C en cloi- 
sonnant la bobine en sections ou galettes montées 
en cascade et enroulées suivant le principe exposé 
plus haut. Un calcul très simple donne alors pour la 
capacité répartie la valeur 


où 7!’ est la fraction de longueur axiale totale 
L occupée par le diélectrique. 


. . I 
En particulier pour + = -; 


» 


Si l'est assez grand devant e la capacité C est donc 


plus petite que dans l’enroulement par couches. 


Pour une gorge carrée 


valeur comparable à (11). 

Avec le cloisonnement, seules les galettes extrêmes 
prennent une charge électrique importante sous 
l’influence des galettes voisines, tandis que les 
charges intérieures se neutralisent deux à deux. 
La capacité propre de la bobine est donc pour sa 


($) En basse fréquence les pertes sont plus faibles et l’inter- 
valle d'air peut être supprimé. 
(7) La capacité radiale des galettes est négligeable, 


20. 


290 


plus grande part rejetée aux extrémités de l’enrou- 
lement. Cependant si la subdivision en sections n’est 
pas assez poussée, l'existence de capacités le long 
de J’axe risque de créer dans chaque cloisonnement 
des déphasages du courant uniforme qui circule 
dans l’enroulement. 2 

Lorsque le coefficient de self est assez élevé 
(L => mH), la technique moderne des bobinages 
utilise des enroulements dits « massés » dont la 
capacité propre est d'autant plus faible que les 
spires qui voisinent dans le bobinage ont un numéro 
d'ordre très différent. En observant cette règle, 
on a été conduit à la meilleure disposition pos- 
sible qui est celle du nid d’abeilles (honeycomb coil) 
dans laquelle les spires chevauchent l’une sur l’autre 
et sont décalées suivant des éléments d'hélice orientés 
successivement de gauche à droite et de droite à 
gauche. fl en résulte donc un croisement des fils des 
couches voisines et une séparation des spires paral- 
lèles par une distance au moins égale au diamètre 
du fil. 

Ainsi les liaisons haute et moyenne fréquence 
des lampes amplificatrices sont constituées de petits 
nids d’abeilles de trois centimètres environ de dia- 
mètre moyen. Par contre pour les ondes courtes, 
on utilise les enroulements cylindriques à une seule 
couche à spires jointives ou écartées. 

En conclusion de ces remarques nous pouvons 


‘ fi r ne El 
dire que la valeur C — _ D U.E.M. représente l’ordre 
de grandeur de la limite supérieure de la capacité 
répartie si le bobinage n'est pas compact ou s’il est 
réalisé sous forme de nid d’abeilles. 
La condition de bonne résonance s’écrira donc 
+107 20 2). 


DES 


I 
9 
On pourra prendre par exemple 


EVE AMC MATOS (49) 

Applications. — Différentes variables qui sont 
l'encombrement, la fréquence, la self, la puissance, 
le volume du métal, le facteur de surtension inter- 
viennent dans la technique pour fixer les dimensions 
et le mode d’enroulement d’une bobine. 

Nous allons étudier quelques exemples. 

Supposons que la bobine soit insérée dans un cir- 
cuit de capacité l bien déterminée destiné à émettre 
une longueur d’onde À. Toutes corrections faites 
pour tenir compte des fils extérieurs à la bobine, 


la self est fixée par la relation L — el 
(on 


10 On se donne D (à partir de l’encombrement 
par exemple). On a 
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= 2 FH) . & ÉD ER 


Cette expression présente un minimum @ pour 


‘ 


a 1}? 
PNR = —— ; ï 
r2nBw?7rt 
soit 
D) Ja 
VND === — 5 
(a TU) 7°? nb (13) 
Comme 
re n D 18 
LTD el Ro Te = rrro/ "6, 


on en déduit 


2L7rw7r 


re 


16 n6 
LL —onrrvo nb, 
o a 


oLrvr fn 
ADO Ve 
La résistance minima en haute fréquence est double 
de la résistance en courant continu. , 

L étant fixé, il y a donc intérêt à utiliser un fil fin 
(2r petit) enroulé sur une bobine de grand diamètre. 
moyen D. Aa 


CE Ro == 


(44) 


De plus le facteur de forme est Ba: Orona 
œ 


el 


CE | B 


r?2D) 
= ——) 


l 


et physiquement il est logique de prendre pour valeur. | 


e 

F. 
Dans ces conditions le facteur de forme est sensi- 
blement indépendant de la forme de la bobine ‘et 


approchée de n, rapport des énergies la valeur 


, x A 5 Q Q RC A “ä 
égal à C'est bien ce qui est vérifié par l’expé- 


rience, du moins pour des gorges à section rectan- 
gulaire (8). ; "+ 
On en déduit la résistance optima 
2Lrvo 


RE 


(15) 
valeur remarquable indépendante de la conducti- 
bilité, donc de la nature du métal. : 

Pour l’optima, l'égalité (13) donne le nombre de 
brins par fil g 


(16) 
et la section 
Diese 


S = rNr= ETS 
r 0) nL 


(17) 
Remarquons aussi que la condition (4) d’unifor- 


2 2 
= TU RS 


Ë 2R V pe 
(en première approximation). Donc le facteur de forme 


2R 


(5) Avec une gorge circulaire : 1=— 
D 


ET ANT US LM 1 
vaudrait —{—) qui est voisin de —. 


ce à 


# 
ui 
ñ 


RÉSISTANCE EN HAUTE FRÉQUENCE 


; mité du champ qui exige que la quantité 


I a D 
Nr\ 4nZ 


! 


k soit très petite donne avec (6) : 


(U) 


7 AW Y7? 
nZ nee AURTR 
5 2 
(4 D? PnE Æ ; 


* Il faut donc pour que la condition (4) soit réalisée, 
£ ane l’on ait 


I 
SVsror TOY 


; | Cette. inégalité exprime que 
1 


| À 


de 


7 - peau e —(2707) * est supérieure au rayon du brin. 
Autrement dit la condition nécessaire et suffisante 
d’uniformité du champ pour la résistance optima est 
# celle de pénétration complète du courant dans chaque 
” brin. 

Cette conclusion est en accord avec notre hypo- 
thèse initiale. 


l'épaisseur 


_ Application numérique. Considérons une 
+ bobine faisant partie d’un circuit destiné à émettre 
lune longueur d'onde À — 100 m, la capacité T totale 
étant 7.10-5uF. On en déduit © —6 7.10 et 


101 
= = - U.E. 
£ 7,307 ce 


“ Supposons l’enroulement en fil de cuivre 


Ja 
_ isolé de diamètre or 


np Un DT, 


- tance optima aura pour valeur 


Le nombre de brins par fil donné par 


À D? VA % 
N= :=—— or 
Z TU y 7 nL 


I 
A D 4 


S— Nr =10,08 mm? 


 L'épaisseur de peau est 


ce qui correspond à une valeur courante. La résis- 
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Cependant l'expérience et la théorie montrent 
qu’au-dessous de room voire même 200 m, le bobi- 
nage à fil divisé est plus résistant qu'avec un fil plein 
de même diamètre. Cela tient à ce que les pertes 
dans l’isolant qui entoure chaque brin deviennent 
trop importantes pour des fréquences supérieures 
à 2 000 ou 3 oookc. 

Prenons un autre exemple 


IX —=1300 M, DNS HE TO LURE 
D'où 
—2 6 
DRE H = UE, M: 
67? ( 
I A : A 
Avec 2r — Far 0n déduit par les formules précé- 
dentes 
% l 
R—= 2,65, IN (pour £ = = ;) 
et 


$ — 0,06 mm?°. 


La condition d’uniformité r <e est évidemment 
satisfaite, 

J'ai dressé ci-dessous un tableau qui pour une 
gamme de longueur d’onde variant de 200 m 


à 20 oo0o m. donne les valeurs de & de N et de s 
k 3 (Tableau I) (°). Les résul- 


| te 2, I, 3? # 
tats obtenus par ce calcul sont en excellent accord 
avec ceux qui sont donnés par les constructeurs 
pour des enroulements du type « soigné » réalisés 
par exemple en nid d’abeilles (self-mignonnette) et 
utilisés universellement à l’heure actuelle. 

Au dessus de 300 m. les valeurs de & obtenues 
avec du fil divisé sont toujours inférieures à celles 
correspondant au même bobinage réalisé en fil 
plein. Remarquons enfin qu’un accroissement du 
diamètre D (lorsque cela est possible) produit tou- 
jours une diminution de &. 


pour 


Remarque. — Pour diminuer la capacité répartie 
il y a intérêt à ce que le volume occupé par le 
métal r?r?nND soit petit par rapport au volume de 
la gorge qui le contient rx D® c’est-à-dire à ce que 
le rapport 
rx r2Nn 
a D? 


« 
OS 


soit petit. A l’optimum on a 


1: D œ 
AS ue TU) 13 F 
Donc 
Ô — (18 
5 oyr A ang ) 
Avec 
el de x? 1) re 
TR Re et 0 LED 


(*) Les valeurs de T et de Z indiquées dans le tableau 


correspondent à des circuits normaux de réception. 


992 


on obtient < I D 
O—=———— —: 
rwyrD e 


En supposant - —{4 (ordre de grandeur assez 


fréquent) la valeur de à s’écrit finalement 
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Cette quantité devant être inférieure à l'unité 
(wy étant donné) il en résulte une limite inférieure 
du produit rD au-dessous de laquelle le bobinage 


devient compact et par conséquent de forte capa- 


cité répartie car le coefficient k de la formule C =kD 


A È Dee , 
D est alors supérieur à —- 
ZWY 7" D) v? 

= | 
| 
TABLEAU Î. = | 
l l u LEA LA ET 
— = 4, — —=,2. — = l ==: ne Es | 
D D D D:2: D 4 DST | 
Ro: ———— — " —— Te RS LE En La 
2m). FmuF) ZL(mH), D=#em. N,. s(mm?). N. s(mm),. N. s (mm”). N. s(mm°). N. s (mm) N. s (mr 
| 
LOOM AT 10 - 10,090 TD K1% 0,2 29 O,14 20 0,1 14 0,07 10 0,05 7 0, 
200  O,11 o,11 2,4 2H ON 1022010 CON 9,5 0,07 6,7 0,05 472810 d 
300 10310 ONTOTAUES 06 31. 0,24 29 SONT FD DAT OLA 11 0,09 5500700 5,505 
600 0,23 0,438 3,49 38 0,30 29} 0,21 19 0,19 TI FO CON TT 9,2 0,07 65708 

000 M0: 0,934. 4,62 43 0,34 30 0,24 D D DOS 1, ÉNEETo E 10,7 0,08 7020 
3000 0,64 HS 201070 62 0,46 14 0,34 31 0290 ONE FO MAO 11712070 
6000 1 10,13 7,9 80 0,63 6 0.44 40 0,32 28 0,22 20 0,16 14 2 epEn 
TODODMMETE 0 18,08 8,85 99 0,77 TO 0,04 49;5 "0,39 35 0,27 20,71 019 17,0 MONT 

D = 6 em. 

13000 2 31,01 6,63 208 1,6 1495: 1,10 104 0,80 74 0,98 02 0,40 J'AI 

20 000 2:14 16,66 T0 228001.%0 1001, 24 11/ 0,80 80 0,62 5 0,41 40 0,3 
É t 


Comme la valeur maxima de r est fixée par la 
condition de pénétration on en déduit donc une 
limite inférieure de D. Et les considérations d’en- 
combrement, la nécessité d’éviter l’étalement du 
champ magnétique extérieur de la bobine, la possi- 
bilité de réaliser de fortes amplifications ont conduit 
la technique à l’adoption de petits diamètres moyens 
malgré l’augmentation corrélative de la résistance 
optima. 


Exemples. — Sur l’onde de 20000 m avec D — 6 em 
oi obtient 


Ô — 0,88. 


L'espacement des spires est tout juste suffisant. 
Sur 15 000 m avec D — 6 cm on a à — 0,67, ce qui 
permet un bobinage à faible capacité répartie. 
D'ailleurs au-dessous de 15 000 m on peut employer 
un diamètre plus faible. Par exemple pour 


À = 10000 m, Di rcmE 6 — 0:07 


Pour des fréquences croissantes d est évidemment 
de plus en plus faible : sur l’onde de 300 m avec 
D°—="{, cm ou D — 2 cmona 


0 — 0,02 ( Ô — 1 
Ô — 0,02 Ju Ô — 0,04, 


les enroulements peuvent être très satisfaisants. 


20 La valeur de la self étant toujours donnée, 
ainsi que le rayon du brin, supposons qu'au lieu 
du diamètre moyen D on se fixe le volume du 
métal V = r’r#nND. Il s’agit de déterminer la 


bobine qui utilise le mieux son métal pour une self … 


et un volume donnés. 


Fr 


Transportons la valeur D) — ue 


T7 *n1 Æ 


mule générale de la résistance (5). On obtient 


ré ri0nfpyw2 V3 75 | 


a? 


L 1 
ae rs 


dont le minimum a lieu pour - 


ne # s É | 2 £ 
ee PEYEETS _ C4 


La résistance optima est alors . 


RE 
( (AO Æ) 3 
e 


D'autre part È 
n D L 


Ro= —= = 
yNr? 


Gyr2nN N 
qui en tenant compte de Nn à l'optima donne 


Ro —= PER NT ET 10 
"14 C 
HE AO) 


GE à Ro. 


Donc 


& est donc plus voisin de la valeur ohmique R, que 


- dans la for- 


à dans le cas DRAEUEnE où l’on s était fixé le diamètre 
moyen. È 
En prenant D comme variable au lieu de Nn on 


la résistance R; passe 


É. aboutit au même résultat : 
- par un minimum & lorsqu'elle est égale aux { 1 ; de la 


: 


a 


: résistance ohmique. Très sensiblement on a 
4 1 
De 5 u) 7 nt 

R= LT: le 
1e F “Gr) (#2) 


1 
% ñ . 
n’est plus approxima- 


(20) 


Er. 
= 


Le facteur de forme { 
4 pa 

tivement indépendant dé la forme de la bobine 
comme dans le cas précédent. En remplaçant CA 


7 n par leurs valeurs il s’écrit (£ I 3): 


# D'où 


(21) 


De plus l'influence du diamètre moyen sur la 
1 


2 


à 
à 
É 
; 
». 


rapport - aussi petit que possible. Cependant / ne 


leur de & est moins importante (D ? au lieu de D-'). 
Par contre, la longueur de la bobine joue un rôle 
important. À première vue il y a intérêt à rendre le 


doit pas descendre à une valeur trop faible afin que 


TANT 


la définition de A soit toujours valable. 
L'étude se poursuit comme dans le premier cas. 


. En particulier à l’optima 
# DV 
D Fe 99 
È nee 4 (z Ha) Se Fu 
= et ; 
rs x 
Dy \” 
= 7 . 93) 
4 à NÉ 2 . Lr ni 5) 


Dans l’ensemble, cette condition conduit donc à 
l'adoption d’une forme de bobine moins massive 
que dans le premier cas. 


= 3° Comme dernière application envisageons le 
cas d’un circuit d’accord formé d’une self L et 
- d’une SPAr F en dérivation. A la résonance 
4 
5. 


bon aT— «+ Supposons que ce circuit reçoive 


7. 

une puissance W. A la résonance, son impédance 

S 24 

6 Le et si E est la différence de potentiel 
F 

pau apparaît aux extrémitès de la bobine, on a 


Rr 
L? 0)? 


W= EE 


Pour une puissance donnée, E sera d’autant plus 
-est petit. D'autre part nous avons 


grand que — 


ET 


vu que l’on pouvait poser L — k D. D'où 


] 
Dos 


; D'UN ENROULEMENT A FIL DIVISÉ 293 


Le problème revient donc à chercher le minimum 
de l’expression 


= 


RF k D 72 nBw2r6 
de le mme 0) 


qui a lieu pour 


D la 


Rent 
zur \V 15 


La valeur de 6 est alors 


6—okrwr(/ À. (25) 
x) 

Avec nos DAS conventions sur &,{G,1n, on 

tire 


e 
= kwr—=92Tw— 10), 26 

2kw7 210 (210:) (26) 
Comme dans le premier cas étudié, la valeur optima 
de 0 s'obtient quand la résistance apparente est 
double de la résistance en courant continu. Le facteur 


de forme (/ -L est le même. 
œù 


Û 
Mais la relation TL = k D exige pour sa validité 
un enroulement peu compact. Cette condition 
conduit à des bobinages moins massifs que dans 
le premier cas. 
Prenons l’exemple d’un bobine répondant aux 
caractéristiques suivantes : 


À = 200 m, & = 0,96 w, D =10;cr»;, 
I j 
L'or rm; 2 = — Cm. 
109 
On a 
= or 15 
Ta 0100 10m: 


D'autre part à l’optima 


I ë £ 
0 = 2 Æwr — 2.107 2037.100 — —0,04:107 15, 
200 


RAT 
Tr 6 © 0,9 est excellent, 


Par contre, un bobinage pour grandes ondes 
À — 20 000 m avec 


Le rapport r — 


R = 4,22 0, D = 10 em, 
A I 
L = 16,66 mH, 27 = —cm 
100 
donne 
R& 
Trw 21 L0mE et 0='0;94 T0 27, 

D'où 


& 
Liv? tw? 
soit environ vingt-deux fois plus faible, la différence 


En modifiant l’enroulement de façon que 


(0) Les dimensions de la bobine interviennent par l'inter- 
médiaire du facteur k, 
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de potentiel disponible aux extrémités de l’enrou- 


lement (pour une même puissance) sera Vas = =—4,09 
fois plus grande. 

J'ai dressé un tableau donnant les valeurs de 7 
pour des bobines de diamètre D = 10 cmet D = {cm 
dont la résistance est celle correspondant à l’optimä 
(deux fois la résistance ohmique) et pour des lon- 
gueurs d’onde comprises entre 200 m et 20 000 m 
(Tableau IT). On voit que + est d’autant plus faible 
que la fréquence est plus élevée et que le diamètre 
moyen est plus grand, c’est-à-dire la bobine moins 
compacte. Il faudra donc agir dans ce sens sur le 
mode d’enroulement si l’on veut recueillir la diffé- 
rence de potentiel maxima pour une puissance 
donnée dans le circuit oscillant. 


TABLEAU II. 
L 
Rio R, 

ÀA(m). L?w? L?w? 0. En: ch 
200 0,89: 1071502192 10712 EN0,04. 1013200027 2,30 
S00 NO OI TOC RO DO TOM NO D TOR SET 10 
C0 MO TS TON NET 00 TOR NO SL MIO 2:50 0,70 

HOOO MO; 0 7 10m TO TOR MO LONTOR PRES 02032 

Glooo on ion 0 770 Tome 0 681010870025 

LNOOOM T7 LOT 20, TORLVEMO 12, 107107 35 
20000 = T0 0254 10m 10 04. 10m 0722 55 


Conclusion. — Notre étude nous a conduit à une 
formule très simple donnant la résistance en haute 
fréquence 


Rr= R Fast ZrovnN * RARU ES 
< x AD A ( 2 ; 


et pour de petits angles 


Rr—= Ro [r+ . (Ground) 1 


Ces formules sont plus maniables que celles de But- 
terworth et d’Austin et donnent des résultats 
comparables. 

Nous avons vu qu’à chaque cas d'utilisation cor- 
respond une résistance optima @ et une forme 
particulière de bobine, toutes deux en parfait accord 
avec les résultats expérimentaux obtenus par la 
technique des bobinages. 
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lisation : encombrement, puissance, self-induction, 
facteur de surtension etc. 
D'une facon générale, pour les solénoïdes à une 
sos L : 
couche, on choisit le rapport 7j Compris entre 0,4 et 1 


si possible voisin de 0,4. Pour une galette plate à 

e I 

une couche, on prend souvent Fe 
4 

selfs plus élevées, bobinées en solénoïdes à plu-. 

sieurs couches, on prend en général les proportions. 


suivantes (11) 


s ; avec l ae 
D 7 ave D — 8 
ou x. 
e : l 
een avec D 
ou 
e : {} : 
Dati avec D — 0,2. 


Remarquons aussi qu'au-dessus de 3 000 ke, le 


fil plein devient préférable au fil divisé, même. 


2 Re, 
Il n’existe donc pas de forme yes de une 
bobine. Cette forme dépend des conditions d'uti- 


Pour des 


«très divisé ». En effet, les pertes dans le diélectrique 


deviennent importantes. 


qui isole chaque brin, 


Îl existe une fréquence critique au-dessus de laquelle 


un fil plein de même section offre moins de résistance 
qu’un fil divisé. Rogowski donne comme valeur 


() 
l’angle d’enroulement). 
Cette formule correspond Re à la réalité. 
C’est la raison pour laquelle la partie haute fré- 


critique de la fréquence ; + 
VER ee 3 3,405 
NE Bor 
_(N nombre de brins, r rayon du brin, ee 60 À la 
température ordinaire et 5 — 3, 54 5 d de 0 étant ñ 


quence des blocs d’accord est constituée par un. & 


petit solénoïde en fil plein. 


(#) Ce sont les dimensions employées dans la construction 


des selfs pour appareils de mesure, ou dans les selfs étalons. 


Manuscrit reçu le 6 avril 1949. 


BIBLIOGRAPHIE. 


[1] LANGEvIN P. — Le champ électrique et le champ magné- 
tique. Cours du Collège de France, 1936. 

[2] SOMMERFELD. Annalen der Physik und Chemie, 
1904, 3, 673. 

[3] MAXWELL J. C. 
tisme, 2. 

[4] CoLomBANI A. —Étude sur les courants de Foucault, 
Journal de Physique, novembre 1948, 9, série VIII. 


— Traité d'électricité et de magné- 


[5] MESNY R. — Radioélectricité générale, 1. 

[6] Bouasse. — Oscillations électriques (Delagrave). 

[7] Buner P. Solénoïdes, écrans et transformateurs 
(Baillère). 

[8] SCHELKUNoOrFF. Electro-Magnétic 
Nostrand Company). 

[9] BuTrerworTx $. — ÆExpérimental Wireless, 3, 1926. 

[10] AusTIN B, B, -— The Wireless Engineer,11 janvier 1934. 


Waves (Ee Van 


SÉRIE VIII, TOME X, OCTOBRE 1949. 


E= TENTATIVE D'EXPRESSION DU PHÉNOMÈNE D'HYSTÉRÉSIS DE TRACTION, 
“% POUR UN POLYMÈRE ÉLASTIQUE CRISTALLISABLE (!) 
; Par Jacques CLOUAIRE. 
Ingénieur A.M., Ingénieur Docteur, Ingénieur du Caoutchouc I. F. C. 
Sommaire. — L’étirement d'un polymère élastique cristallisable (tel que le caoutchouc naturel vulca- 
+ nisé, type « pure gomme ») dans les conditions habituelles des dynamomètres commerciaux, présente 
“4 les caractères suivants : 


/ # . DE £ L: A La 
1° Irréversibilité thermodynamique des phénomènes. 
20 Liaison héréditaire entre la valeur de la tension à un instant donné et les déformations antérieu- 


rement imposées. 


On exprime la tension résultant d’un allongement donné, à un instant donné, en tenant compte de 
l’élasticité de la phase cristallisée, de celle de la phase demeurée amorphe, et de la vitesse de cristalli- 
sation à cet instant; un terme additionnel tient compte de ce que l’on appelle le frottement interne. 


I. Caractères expérimentaux de l’hystérésis. 
— Des essais expérimentaux (1) portant sur l’étude 
de l’hystérésis développée par des échantillons 
- de caoutchouc naturel vulcanisé, soumis à des 
 tractions et rétractions alternées, l’allongement 
pouvant aller jusqu’au voisinage de la rupture, 
_ ont permis de dégager les caractères principaux 

- de ce phénomène, et, à leur Jumière, de tenter 
l'interprétation qui fait l’objet de cet article. 


Les conditions opératoires furent les suivantes : 


1° Utilisation d’éprouvettes Schœpper standard; 


- 20 Tractions et rétractions effectuées sur dyna- 
- momètre Lhomme et Argy; à une vitesse d'extension 
_ de 30 cm: min; 


30 Tractions et rétractions se répétant sans inter- 
ruption susceptible de provoquer des effets de 
reprise élastique; 


4° Cycles décrits pour toutes les amplitudes pos- 
 sibles, exprimées selon un multiple entier de 100 
pour 100 d'’allongement relatif, pour des allon- 
 gements compris entre o et 700 pour 100 (voir 
tableau I); 


50 Nombre de cycles le plus élevé possible, ne 

* dépassant pas cependant 2000, mais plus limité 

* en général (50 à 200), à cause de la rupture pré- 
_ maturée des échantillons. 


(4) Ce Mémoire développe les principales conclusions 
* d'ordre théorique auxquelles ont conduit les recherches 
eflectuées par l’auteur à l’Institut Français du Caoutchouc. 


: 


Les résultats obtenus montrent que la nature même du phénomène de cristallisation par étirement 
(tendance vers un état de cristallisation partielle, dépendant de l'allongement imposé, et de toutes les 
déformations précédemment subies) explique le phénomène d’hystérésis de traction particulier à ces corps. 


“ 


60 Conditions de température ambiante et d’état 
hygrométrique controlée à 209 + 10C, et 6o pour 100 
10 pour 100 (ces valeurs n’entraînant pas de 
variations décelables des quantités mesurées) (?). 


TABLEAU I 


Allongement maximum 
0% | 100% | 200% | 800% | 400% | 500% | 600% |700% 
ON] 9.0] 0-10 0.200 0: 800 0:40 | 9. SN ]0. 600 | 0. r0ù 


Ù 
10e 100.100 100-200 100-300 100-400 ne 500 00.600 100.700 


200% | 
200.200 200.300 200.400 RUE 500 
300% oo 
300-300 (800.400 300.500 (300.600 1500. 700 


400 % 


MINIMUM 


400.400 400.500 400.600 400.700 


500 % 
(500.500, |500.600 [500 700 300% 


Le | 
600 % 
A) 600.600 |600. 0 | ; 50% 
700% 
700.700 100% 
Essai Lis 


Allongement 


Les observations peuvent se résumer aux principaux 
faits suivants : 


19 L'’hystérésis, mesurée comme l’aire d’un cycle 
(exprimée en unités d'énergie) rapportée à l'énergie 
fournie par la machine dynamométrique à l’éprou- 
vette, lors de la traction, croît avec l’amplitude 
des cycles. 


(2) L'essentiel des résultats d'essais obtenus a été publié 
dans la Revue Générale du Caoutchouc à laquelle on pourra 
se reporter pour plus amples détails expérimentaux. 

J. CLOUAIRE. — Rev. Gén. Caout., 1949, 26, 85-90. 
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20 Pour une amplitude donnée, (Healsee dans le 
domaine des allongements de façon définie, l’hys- 
térésis diminue lorsque-le nombre de cycles augmente, 


30 En général, les cycles successifs ne se super- 


posent pas les uns aux autres, mais se déplacent 


vers les tensions plus basses. 

Bien qu'on n’ait pu mettre en évidence une limite 
à cette migration, son importance diminue lorsque 
le nombre de cycles augmente. 


4° À amplitude égale, l’hystérésis est sensiblement 
indépendante de la localisation dans le domaine 
des allongements alors que la migration augmente 
lorsqu'on se déplace vers les allongements plus grands. 


50 La migration est d’autant plus accusée que 
l'amplitude est grande. 


60 Le repos prolongé entre les cycles successifs 
entraîne le retour lent vers les valeurs déjà observées. 
Ainsi, pour des cycles décrits entre o et 700 pour 100, 
le premier et le deuxième étant décrits sans inter- 
ruption, et le troisième après un repos de durée 
variable de l’éprouvette, en l’absence de toute tension 
imposée, au bout de 8 jours, ce troisième cycle est 
superposé au deuxième. Au bout de 30 semaines, 
ce troisième cycle est compris entre le deuxième et 
le premier. Il semble toutefois que le retour aux condi- 
tions exactes de l’origine soit très lent, sinon impos- 
sible. 


7° Il semble a priori que l’on pourrait délimiter 
la zone d’allongements, de o à la rupture, en deux 
régions : 

a. au-dessous d’un allongement sensiblement égal 


à 200 pour 100, l’hystérésis est faible, et la migration 
négligeable. 


b. au-dessus de cet allongement, la migration 
s’accuse suivant les principes mentionnés ci-dessus, 
et l’hystérésis augmente considérablement. 


On rapprochera de ces faits les observations 
suivantes : 


a. le volume d’un mélange de caoutchouc du type 
gomme pure, vulcanisé, reste sensiblement constant 
jusqu’à un allongement de 200 pour 100 environ, 
et varie ensuite (1), (2), (3), (4). 


b. les diagrammes de rayons X montrent que la 
cristallisation devient notable, pour des allongements 
supérieurs à 200 pour 100 (5), (6). 


c. des diagrammes de rayons X ont montré que, 
pour un allongement donné, la cristallisation se 
poursuit dans le temps, lors même que la tension 
diminue (7). 


d. les éprouvettes tendues dans les conditions 
définies précédemment exhibent une augmentation 
de température notable (de 5 à ro°C) pour des 
allongements supérieurs à 200 pour 100 envi- 


ron (1), (8). 
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+ 3. Tentative 
mènes. — A. CAS D'UN POLYMÈRE NE CRISTALLISANT 


_e. les théories statistiques de l'élasticité, condui- 
sant à des équations comme celle de F. T. Wall, . 


par exemple, ne sont sensiblement vérifiées par À 


l'expérience que pour des allongements au plus 
égaux à 200 pour 100 (9). 


2. Conditions d’ élasticité des hauts poly- 
mères cristallisables. --— Loin de considérer, 


<-œH 


dx 


comme on l’a fait parfois, l’hystérésis comme une 


irrégularité gênante, nous la tiendrons comme une. 


des manifestations fondamentales de la nature de 
l’élasticité des hauts polymères cristallisables, aux. 


grands allongements. A la lumière des Re 
citées, on peut dégager les deux caractères fonda 
mentaux suivants : 


10 Les conditions thermodynamiques de l’élas- 


ticité s’écartent d'autant plus de la réversibilité 
que l'allongement est grand. 


20 La valeur de la tension obtenue, pour un. 
allongement donné, dépend de toute l’histoire. des | 
déformations précédemment imposées. 


tères seront que : 


_19 On ne pourra en général utiliser la notion de : 


% 


< 


Les conséquences immédiates de ces Re carac- È 


réversibilité dans le traitement hernie 


des phénomènes. 


29 On ne pourra utiliser des équations algébriques 
classiques dans la traduction mathématique des 
phénomènes. Si y est la fonction étudiée, et æ, ds 

., 4, les variables indépendantes reconnues | 
comme agissantes on ne pourra former une équation 


La, 


ME DIU 00 


Zn) 


telle que pour des valeurs de % Lys Lys + 
nées, y se trouve déterminée. 


d'interprétation des phéno- 


PAS PAR. TRACTION. — Nous admettrons en général 
qu'il existe une certaine hystérésis, attribuable à 
ce qu'on appelle « frottement interne » (#). Dans ce 
cas, la réversibilité thermodynamique ne pourra 
s'appliquer, et nous aurons seulement d’après a 
premier principe : 


FA = dU, —dQ, (l'indice a étant utile par la suite), 
F étant la force extérieure appliquée, l la TL 


de l° éprouvette (nous négligeons ici le travail de la 
pression hydrostatique extérieure). 


En raisonnant sur l’unité de volume, nous aurons 


, Ln AOD- 


] 


(5) Sans entrer dans le détail des causes structurelles de ce 


phénomène, nous appellerons frottement interne l'effet dont 
résulte une hystérésis cyclique (les cycles successifs se super- 


posent à eux-mêmes), en l'absence de toute cristallisation. 


C’est ce qu’on observe expérimentalement, lorsque les allon- 
gements sont inférieurs à 200 pour r00. 


> ee 


(43: A5) 


3 5 étant alors la tension par unité de section initiale 
M (avant toute déformation) et 2 l'allongement relatif 7 


Nous considérons que ce que nous appelons Us 
Ron interne n’entraîne aucune variation de l’éner- 
“gie interne pour un processus isotherme, c’est-à-dire 
qu'il s’agit d’un phénomène parfaitement identi- 
fiable au frottement, une partie de l'énergie fournie 
par la machine de traction se transformant en 
chaleur, Soit Q, ia Gusnuse de chaleur ainsi déve- 
_loppée. Il vient : ; 
É Éd 2 die 00: 
- | ut 
24 | 
Lorsque Q, — o il n’y à pas de frottement tentes 
et on a affaire à l’élasticité parfaite. Nous pourrons 
alors considérer que la réversibilité thermodyna- 
._mique s’applique. Si l’essai est conduit assez lente- 
- ment, l’isothermicité sera respectée : c’est ce qu'on 
observe pour les essais conduits aux vitesses clas- 
- siques des dynamomètres, 10, 20, 4o Cm : mn, pour 
des allongements inférieurs à 200 pour 100. 
Nous pourrons alors remplacer les dérivées totales 
par les dérivées partielles, à température constante 


PRE OU, 002 CION , Re 
SE L( OX —. OX ) | (Fe 7+ )» (2.3 À) 


ou encore x : 
AONACECRES 


Lorsque (5) — 0, la tension est déterminée, au 
T 
=. 
Ô 


OX - 
(facteur d’hystérésis), par le 


28 
ox 
utilisant l’équation de Boltzmann, reliant l’entropie 
à la probabilité thermodynamique d’un état, conduit 
aux expressions de la tension dérivées des théories 
statistiques déjà mentionnées. 

Cependant, pour des allongements assez élevés 
(supérieurs à 100 pour 100) et en l’absence de toute 


r la transformation de ce terme, en 


cristallisation, il devient clair que sous l'effet du 
- champ de forces internes créé par l’application d’une 

- tension axiale, l'orientation des éléments du réseau 
4 qui forme l’élastomère ne se fait plus au hasard, 
. mais tend à se produire dans une direction préfé- 
_rentielle. 

L'arrangement au hasard des éléments d’une 
5e lorsque seule la distance entre extrémités 
se trouve imposée, indépendamment même des 
- empêchements stériques, est de moins en moins 
- possible. De plus, dès qu’une portion de chaîne se 
: trouve orientée suivant l'axe de traction, toute 
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extension ultérieure ne pourra se traduire à son 


égard que par une résistance de nature géométrique 
(déformation des angles de valence, augmentation 
des distances entre atomes de carbone, ou simple- 
ment déplacement relatif des éléments déjà orientés 


les uns par rapport aux autres, pour des chaînes 


suffisamment voisines). Tout cela se traduira par 
une augmentation de l’énergie interne, au détriment 
du travail mécanique fourni par la machine de 
traction. 

Les limites extrêmes correspondront donc : 


a. au cas de l’élasticité statistique pure 


ESRI JS, jrs 

Ta TS): (4.3 À) 
b. au cas de l’élasticité classique : 

PL Ua: L € 

= ( + ), (5.3 A) 


l'indice { indiquant que U ou $, joue exclusivement. 
Pour un état d’allongement donné, nous pourrons 
dire que 7, contribuera à la valeur de la tension, 
d'autant plus que l’allongement sera faible, et que +, 
y contribuera d'autant plus que cet allongement 
sera élevé. 

Nous caractériserons ces contributions relatives 
par l’emploi d’un certain coefficient « représentant, 
si l’on veut, la fraction de section d’une éprouvette, 


qui résiste à cause de ty, (1 —) représentant la frac 


tion qui résiste à cause de Tas De sorte que nous 
pourrons écrire : ; 


ra OU { dS2: 90: Q., 
ef) 6 TS },= 5), (RE 


équation dont les termes s’identifient à ceux de 
de (3.3 A) suivant : 


 (OUas Ua 

Où Jr \ dù 

DEN OS: 
or =r(), 


On aura donc un moyen de traiter de l’élasticité 
des corps hautement élastiques, demeurant amorphes, 


par l’étude séparée des conditions aux limites, qui 
fournira le terme en S,, et le terme U,. Le terme 


en Q, se déterminera par la nature du frottement 
interne (assimilable avec une bonne approximation 
à un frottement liquide). 

Pour la détermination de «, on pourra, pour 
simplifier, supposer des chaînes formées de chaînons 
faisant entre eux un angle constant, et doués d’une 
rotation parfaitement libre. La connaissance de la 
fonction de répartition de la longueur de ces chaînes 
permettra de déterminer, pour un allongement 


donné, la proportion de ces chaînes ayant atteint 


leur extension maximum, par le procédé cinétique, 
et qui, à partir de cet allongement, résisteront avec 


“ 
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augmentation d'énergie interne. On pourra prendre 
cette proportion comme valeur de x. 


B. CAS D'UN POLYMÈRE CRISTALLISANT PAR TRAC- 
TioN. — Nous admettrons que la traction fait naître 


dans l’élastomère une cristallisation, d'autant plus : 


importante que l'allongement est élevé, mais que 
celle-ci ne peut atteindre une valeur d’équilibre 
instantanément, et se trouve affectée par tout le 
processus de déformations imposé avant d'atteindre 
l’allongement considéré; de sorte qu’à cet allon- 
gement correspondra une infinité de valeurs du taux 
de cristallisation, puisqu'il existe une infinité de 
possibilités quant à la loi de déformation. Nous 
dirons donc que le taux de cristallisation est doué 
d’hérédité. 
Soit encore la tension : 
He dUrd0s 
der EUX 


Nous appelons U, et Q, les valeurs de l’énergie 
interne et de la quantité de chaleur échangée avec 
l'extérieur, dans le cas où le polymère reste amorphe 
(cas traité en 3-A). 

Nous appelons U. et Q. les valeurs des mêmes 
quantités, dans le cas où le polymère est entièrement 
cristallisé. Enfin, U, et Q, rendent compte des 
variations des mêmes quantités, lorsque la cristal- 
lisation se produit. 

Nous aurons alors : 


dU, S O, 
JO PAS 0, 


dU, < 
do, >= 0 


En effet, la cristallisation à un allongement donné 
tendant vers un équilibre, cette tendance entraînera 
une diminution d'énergie interne (ce qui entraînera 
une diminution de la tension), alors que l’extension 
de la fraction cristallisée se fera suivant les lois 
classiques de l’élasticité, avec Eia à de 
l'énergie interne. 

Soit 5 un coefficient défini selon les mêmes condi- 
tions que «, caractérisant cette fois la fraction de 
section de l’éprouvette qui peut être considérée 
comme résistant suivant les lois de l’élasticité 
‘ cristalline, (1 —6) caractérisant la fraction qui résiste 
suivant les lois de l’élasticité amorphe. 

La contribution de U, et Q, sera donc propor- 
tionnelle à (1— £); celle de U. et Q,, proportionnelle 
à 6; quant à U, et Q, à un instant donné, l’accroisse- 
ment de leur contribution sera proportionnel à l’ac- 
croissement de 5, de sorte que la contribution totale 
à un instant donné sera : 


es dU, d@ ae 
. CRUE ) 
she df dus d@ d} 
IG REUT FI =) 


du, > 0 
dO,> 0. 


ou encore 


Meoub) 
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Nous aurons sn finalement : re 
AD (TS s =) 
(Te — de) dÀ 
a 
Nous exprimerons — et ÎQa 


dÀ $ 
et Q4 (l'indice # nr. qu'il s’agit de la contri- 
bution totale de U, lorsque (4 — 1) —o ou 1 de Q, 
lorsque « = 0) 

Ici, on ne peut admettre que l'essai est isotherme, 
puisque dans les conditions opératoires, l’augmen- 
tation de température de l’éprouvette est notable. 

Nous aurons donc: % 


“dUz Es {OU OU: d7' 
mr) Gran] 6 
et “] 
dO A 902 JO a7’ € 
dx =] ox ),+ É ee 920 
(ete 
OT <), représente Ja chaleur spécifique de la phase 


amorphe, à allongement constant, soit C qui est 
OUat 

égale à (SA), puisque le travail fourni dans ce cas 

est nul; GC; est en général une fonction de À, D’autre 


part (SE), est égal à T(SE): 


Nous aurons donc pour 7 : 


)Uy 8: 
ohe ea A ne), 
ae Car cu 


d2 
SFR dr à dx 


pe ge 


NN FR UE _(5.3B)% 


Èe 


représentant la part du frottement interne 


ee écrivons Q; pour différencier cette valeur de Q 


valable en fraction amorphe seule). 


L'équation (5.3B) pourra se simplifier; nous 
admettrons par exemple que la fraction cristallisée 
résiste grâce à une augmentation dé U., de sorte 


dO. 

que ee = 0 

10, , 
: _ sera déterminé par les conditions de li essai, ee 


se déterminant par la cinétique de la cristallisation, 
une partie de l'énergie libérée sera dissipée en cha- 
leur, compte tenu des dimensions de l’éprouvette, 
de sa conductibilité thermique dans l’état considéré, 
des conditions de température et de convexion 


TS 
extérieures, 


Le reste contribuera à diminuer la 
_ tension, 


Nous en venons maintenant à la détermination 


. de 6. Nous avons vu que le phénomène de cristalli- 


sation était de caractère héréditaire. Nous dirons 
. donc que 6 dépend de l’histoire des déformations 
précédemment imposées, autrement dit, de toutes 
. les valeurs de l'allongement précédemment rencon- 


“trées (4). 


Appelons { la valeur du temps au moment où 


._ l’on désire connaître la valeur de B, et r la valeur 
‘ courante du temps. 


La valeur de 5 dépend de l'allongement considéré, 


| et de toutes les valeurs déjà prises par cet allon- 
_ gement. À 


: _ allongements Âs A5 


l ( l 
B= fr) + Fa, Xe, hs, mn), 


== Eee 20 


les signes jee ke ..., représentant l'influence des 


.., au temps f{, supposant que 
44 on remonte cat dans le temps. 
C'est-à-dire que l’imposition au temps r et pen- 


dant une durée dr, d’un allongement 2,, entraînera 
. une certaine valeur de B dont il restera au temps { 
. un certain résidu, que nous définirons en valeur 


\ 


2 


spécifique par o(é, r,) dr. 
Nous supposerons que le résidu total est, de plus, 
une certaine fonction de À (r),c’est-à-dire, par exemple, 


que le résidu des actions antérieures est d'autant 


_ plus grand que les allongements imposés furent 
- grands. Au temps {, le résidu des actions antérieures 
sera la somme des résidus des actions élémentaires. 


.l 
R= | F[o(t, r), (r)]dr (6.3B) 


Supposant que toute action héréditaire antérieure 


à un certain instant o pris pour origine est négli- 
_geable, nous aurons pour £ 


É; 


B= fn + f. Flett: Fr) X()l'dr, 


» A(r) représente la loi de traction imposée; dans le 


cas des dynamomètres 


classiques, l'allongement 


est proportionnel au temps de sorte que l’on aura 


fe 


TRES te 
= in 2(hm nm) x 0) | 
lon 


À (r ) 


avec 


on allongement maximum imposé, 
À, allongement minimum imposé, 


 O période de la sollicitation. 


Il sera plus aisé de remplacer cette expression 


(4) Sur les phénomènes d'hérédité, voir V. VorrERRA [10], 


Peu: [12]: 


TENTATIVE D'EXPRESSION DU PHÉNOMÈNE D'HYSTÉRÉSIS DE TRACTION 
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par une série de Fourier, c’est-à-dire : 


— Àm 


Mr 


COSI07 


A 
0 < 
ARRET ESS A7. 


Finalement 5 se composera d’une fonction /f(2,), 
fonction algébrique simple de l’allongement, corres- 
pondant à la cristallistion apparaissant lors d’un 
allongement instantané, et d’un terme complé- 
mentaire, de nature héréditaire. 

Si le taux de cristallisation, à une température 
donnée, tend vers une limite uniquement fonction 
de l’allongement, l’équation intégrale 


l 
[tee r) acmar 
0 
doit tendre vers une limite unique pour { — « quelle 
que soit l’histoire des déformations, c’est-à-dire 
quelle que soit À(r). 
Il semble toutefois que la résolution générale de 
l'équation soit compliquée. 
On obtiendra une première simplification en 
admettant que le résidu d’action au temps r est 
proportionnel à À(r). C'est-à-dire 


l 
R(E= f o(6, 72 
û 


C’est une équation de Volterra linéaire de pre- 
mière espèce. Cependant, le problème classique de 
résolution de ces équations est de trouver À(r), 
connaissant o(£, r) et R(#). 

Ici, nous cherchons à déterminer R(#), connais- 
Sant À(r), mais ignorant la forme de o (£, r). 

Nous pouvons, a priori, fixer quelques-uns des - 
caractères auxquels doit répondre le coefficient © (£, r). 

On démontre que, lorsqu'il existe un cycle d’hys- 
térésis unique (hystérésis magnétique, élasticité 
visqueuse. pure), © est une fonction de (fr); c’est 
ce qu’on appelle l’hérédité cyclique. Ici, nous avons 
déformation et migration des cycles, de sorte que 
les importantes conclusions dégagées pour l’hysté- 
résis cyclique ne pourront s'appliquer. 

Mais, si les phénomènes tendent vers une limite 
unique, © se réduira à une fonction de ({-r) lorsque 
l’on aura { — ©. De plus, la migration des cycles 
montre que 5 augmente constamment, pour un 
allongement donné, c’est-à-dire que © sera néces- 
sairement positif, d’une part, et sera fonction crois- 
sante de &. Enfin la disparition de l’action hérédi- 
taire dans le temps (mise en évidence par l'étude 
de l’influence du repos entre chaque cycle) montre 
que » doit être fonction décroissante de (fr). 


r) dr. (8.3B) 
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4, Conclusion. — [L'expérience montre que, 
lorsque l’on étire un polymère élastique cristallisable 
dans les conditions correspondant à celles des dyna- 
momètres classiques, les caractères principaux de 
l’élasticité mise ainsi en évidence sont : 


19 Irréversibilité thermodynamique des phéno- 
mènes ; 


29 [Liaison héréditaire entre la valeur de la tension 
à un instant considéré, et les déformations anté- 
rieurement imposées. 


Partant du principe de conservation de l’énergie, 
on a introduit une quantité de chaleur Q,, carac- 
térisant l’hystérésis due à ce qu’on appelle frottement 
interne; un coefficient x caractérisant les contri- 
butions relatives de l’énergie interne et de l’entropie, 
pour l’élasticité de la fraction demeurant amorphe, 
ce coefficient étant de nature algébrique simple 
(c’est-à-dire que, pour une valeur donnée de l’al- 
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longement À ete une seule valeur possible 4 és 21 $ 
un coefficient B caractérisant les contributions 
relatives des fractions cristallisées et amorphes, 1 
ce coefficient étant de nature héréditaire (c’est-à= 
dire que, pour un allongement donné, l’importance 
de la cristallisation dépendra de toute l’histoire 
des déformations antérieurement imposées). 

C’est de la forme du coefficient d’hérédité o(£, PA 

intervenant dans l'expression de 5, que dépend 
l’hystérésis considérable qui se développe au- . 
dessus de l'allongement moyen de 200 pour 100 
(allongement à partir duquel la cristallisation se 
manifeste clairement), ainsi que les différents carac- 
tères liés à la migration et à la déformation des 
cycles successifs. 
. L'action de l’amplitude et de la localisation des 
cycles dans le domaine des allongements est inter- 
prétée en ce que la loi d’allongement imposée 20). 
intervient également dans l'expression de p.. 


$ £ 
Manuscrit reçu en mars 1949. 
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SUR LA CORRECTION RELATIVISTE DANS QUELQUES PROBLÈMES DE DIFFUSION 


Par MiRceA DRÀÂGANU. 
Cluj (Roumanie). 


= 
Sommaire. — Dans les lignes qui suivent nous faisons quelques remarques sur l'extension relativiste 
des quelques équations qui servent pour étudier, dans les différentes théories, les phénomènes de 

diffusion. 


_ 1. Il me semble utile pour plus de clarté, de On emploie souvent la notation plus commode 
rappeler d’abord quelques définitions fondamen- {3} 2, 0, 0,, 6, t)dæ dy ds déx des de, 

-tales (1). S’il n’y a ni création, ni anéantissement = n(x,y,2,0 fa, av 0 vat)dædydzderde, de, (3) 
de particules, pendant ou à cause du mouvement, À ; Re 

bot si N est leur nombre total, dont le centre se ©U encore la notation symbolique équivalente 
trouve à un moment donné { effectivement à l’in- pee 4) Fo res ne 1) me NS éà 
_térieur du volume fini V, on peut définir une densité Fe RO 4 
_corpusculaire (moléculaire dans le cas où les parti- pour indiquer la probabilité de trouver au moment { 
D sont des molécules). Une telle une particule de l” CDR Dee dont la vitesse 
densi définie alors par la valeur limite du 

rapport entre le nombre des particules et la grandeur 
du volume, quand le volume devient infiniment 
petit, si cette limite existe. On est donc conduit à 


= 
serait comprise entre D et v + dë à l’intérieur du 
volume d V. 


2, Soit maintenant à un moment donné un 


écrire nombre AN très petit, bien déterminé et fini des 
: pe (1) particules. Ce nombre étant invariant par rapport 
À Ar >0AF au changement de la vitesse des particules dans un 


système de référence déterminé, il s'ensuit qu’il est 

invariant par rapport à une transformation du 
groupe de Lorentz (?). Cette proposition peut être 
démontrée aussi directement. En effet, comme AN 
représente le nombre des particules contenues dans 
l’élément de volume À V — dx dy dz, on peut écrire, 
d’après la définition (1) 


Cette définition implique, comme on l'a déjà 
remarqué, certaines difficultés, car en général, ni la 
valeur de n, ni celle de V ne peuvent s "exprimer 
toujours par des nombres entiers. Tout de même, 
dans l'interprétation probabiliste, il est permis de 
définir le nombre AN des particules contenues 


-dans l’élément de volume AV par le produit nA V, AN > 
en d’autres mots de définir une probabilité de MUR (5) 
présence d’une particule (son centre ou bien elle- 

même, si elle est physiquement assimilable à un Dans un autre système de coordonnées 


EP: matériel) à l’intérieur du volume A V. S'(0',x, y, 7, ), 2 nes rapport à S(O, x, y, z, 1) 


D'autre part, une particule considérée quelconque de la vitesse constante v, on peut définir de même 
prise au hasard, possède par hypothèse aussi une une densité 
vitesse de translation, c’est-à-dire une énergie A (6) 
cinétique. On définit alors en mécanique statistique A7'50 AV 


a probabilité que la valeur de la vitesse d’une A1 = 4» dy! dz’ étant l'élément de volume tridi- 
molécule ait sa valeur comprise entre v (de Com-  yensionnel. Soit w, la composante de la vitesse des 
“posantes v,, v,, v.) et v+dv (composantes particules AN par rapport à l’axe de x dans le 
Ds + du;, v, +du,, v. + dv.) au moment Ê par Système de référence $S. On peut écrire alors la 


l'expression relation (). 

À : | > J À I (TS 

JS, Vs 25 Vas Vys V2 t) de de; des ee v t) du, (2) / (: FA . }r. (7) 
Sr $ L p? LS 

4 | o 

M ction (dy 20e b,,-0., 1) étant celle qu’on (ee 


3, ent la loi de répartition des vitesses. ne 
_appelle SET P (2) M. von LAUE, Die Relativilätstheorie, 1, Braunschweig, 


1921, P. 107. 


() Voir sur ce sujet J. H. Jeans, Théorie dynamique des 
gaz, Paris, 1925, traduit par A. CLERC, D. 15-17. 


(#) J. BECQUEREL, Le principe de relativilé el la théorie de. 
la gravilation, Paris 1922, p. 94; M. von LAUE, loc. cit., p. 108. 
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ou encore, en remplaçant les valeurs (5) et (6) 


AN! [ y 0 RDA, 
RS brie 1 4 Q 0) 
0? (: Ge? ) im, AV RCE 


lim — ï = 
Ar rs-0 AV / 
I 


et, d’après les transformations de Lorentz 


dé = (-"%#), (9) 
dt / 2 c° 


l' étant le temps mesuré en S' et { le temps mesuré 
en $. On obtient alors de (8) 


AN” CCS TTAUAN 


lim Ms er 1M *— 
"Aayrrs.0 AV dé Ars.0 AV 


(10) 
et puisque l'élément de volume quadridimensionnel 
est un invariant par rapport aux transformations 
de Lorentz, AV'df' — AV dt on conclut 

(11) 


AN'= AN. Cr OF -D: 


\ 


Le nombre de particules AN est donc un invariant. 


3. Désignons maintenant par v la vitesse des 
particules AN par rapport à un système de réfé- 
rence S(O, x, y, z, D et soit A V, l’élément de volume 
très petit, mais fini mesuré dans le système propre S, 
concernant ces particules. La grandeur de l’élément 
de volume dans le système $S s’obtiendra par la 
transformation 


san / (£)- 
pe 


Si donc nous définissons par (5) la densité corpus- 
culaire des particules dans le système $S et par 


(12) 


. AN 
To = lim 13 
ÿ Ars o Fe ) 
la même densité dans le système de coordonnées S,;, 
on aura la relation 


É + AN 
nn — AU A7 Us 
r> F>0 
0 4/ LS (2) AV 
I : AN I 
= Re || ) 
/ (5) Ne AV va 2 él ni 
1 | — 1— (— 
C2 c? 
c’est-à-dire 
n = sa) 3 (43) 


On définit encore en relativité restreinte le quadri- 


= pè+ pi + VE, MERDE 


L'élément d'extension en phase invariant par 
rapport au temps (5) s'écrit en mécanique statistique 
relativiste de la manière suivante : 


do = dx dy dz dQx dQ; dQ:, (17) 
ou symboliquement 
S 
de = 4740. (18). 


On posera encore dans la suite pour abréger l'écriture 
Q= Q02+ 074 07: (49) 


D'après ces définitions la probabilité de présence 
au moment { d’une particule en mouvement dans le 
volume à six dimensions (17) sera donc donnée par 


W(x, y, 3, Qx, Or, O2, t) dx dy ds dQ:dQ;dQ:, (20) 


ou bien symboliquement 


w(r, O, :) dy d6, (21). 


W étant une fonction de probabilité généralisée. 
Si Æ représente la fonction généralisée de répar-. 
tition des vitesses (6), on a 


rs W(x;, Ÿ> 3; Q>, Q);; Q £) 


= na, y, 2 t) Fa 7, 2, On On Ont) (O2) 


On peut alors écrire, en tenant compte de (12) 


nr, @,r)araô= (7 0,1) ar, 48, (3) 


car ‘on a 


(?) Voir par exemple J. BECQUEREL, loc. cil., p. 114. 

(5) M. PLANCK, Verh. d. Deulsch. Phys. Ges., 1906, 8, 
D. 136-1413 F. JÜTINER, Ann. d. Phys., 1911, 4e série, 34, 
p. 826-880. 

(9 F. JÜTINER, Ann. der Physik, 4e série, 1911, 34, 
p. 826-882; Ann. der Phys., 4° série, 1911, 85, p. 145-161; 
LOT Physik, 1928, 47, p. 542-566; voir aussi A. BUSEMANN, 
Physik Z., 1932, 88, p. 775-777; P. E. Bônmer, 1928, Æ7, 
p. 542-566; Phys. Z., 1933, 84, p. 454; A. BôHMER et 
A. BUSEMANN, Physik Z., 1934, 35, p. 350. : 


si l’on désigne n, F par W,. L'expression (21) est 
_donc invariante par rapport au groupe de Lorentz. 


- 4. Nous introduisons 


maintenant le quadri- 

vecteur M,, dont les composantes sont 
Mx=W; Qx; M, =W Q), | 95) 
M =W0% Mi=WQu | SEA 


L’équation de continuité de l’hydrodynamique 
_relativiste s’écrira alors de la manière suivante : 

DivM,= 0 (26) 
ou encore 


W Div Qu+ (Qu, Grad Wo) = 0 (27) 
Les opérateurs Div, Grad, appliqués, soit au quadri- 
- vecteur P,, soit au scalaire U, signifient les opéra- 
tions 


OP: De Ps, op, 


% DivPy, = 


en — — jC cs PA 

NE Des D À ou = ic ot, (28) 

nr Le Le To se) (9) 
dy” 03° du 


Tenant compte de (24) et se et en introduisant 
la nouvelle fonction 
enr Qi), Lu 


wo (rs 
Ale p, Je =W, da Q, e), 


| ja formule - (26, 27) Én s’écrire sous la forme 
vectorielle suivante : 


S SRUVS 
d æo (7, P, € . lo mo\r, v, 4 


nn va su + div V = 
el 0 
A c? c? 
où le symbole Div a la signification habituelle du 
calcul vectoriel. La fonction W se transformant 
. d’après la loi (24), elle doit avoir vraiment la signifi- 
cation d’une probabilité de présence, comme l’a 
- bien remarqué E. Durand (?). 


(30) 


OO L) 


5. Supposons les particules soumises au cours 
de leur mouvement à des chocs élastiques. Soit alors 


| po. Q—0Q)a0 
Fe ue Q; Q2 O5 Qx, 
Qy— Qys Q:— 02) dOx dQ, dQ:; 


la probabilité par unité de temps pour qu’une 
particule change sa vitesse Q’, Q, Q: et passe 
per suite d’un choc dans l’élément d’extension 
- en phase dx dy dz dQ; d0Q, dO. 

Désignons par T'(Q:, Q,, Q-) la probabilité totale 
. (9) E. Duran», Journ. Math. pures et appliquées, 1946, 
25, p. 181. 
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(32) 


303 
de déflection à la suite des chocs 
l'(Qzx Qy> Q:) 
= f&(0.0,01,0.—0. 
Py— 0, 0:—0Q:)d0:d0;dQ:, (33) 


puis soit Q (0, Q,, Q-) la valeur finie de la proba- 

bilité totale de capture par unité de temps. 
Tenant compte de (25), (28), (33), (34), la varia- 

tion totale du nombre des particules dans l’élément do 


et l’intervalle temporel df est alors 


> 
(Div M,.) dd dé 
Ÿ > 
— 4Fd0 mr, 


+ av d0 4 W(, o e)Iro)+a(o) 
ou encore 
: ne 
Dix (Qu Wa) = (7, @1)[N(Q)+0(0) 
> > 
+ fr 0,t)%(0,0— Qad. (3) 


Q,t)&(Q Q— 040 


(34) 


On démontre ensuite facilement la formule ($) 


on 
e 
: 


(36) 


D = Qi+ 0$+ 02,1 0 v2+02+ 02, 
La relation intégro-différentielle (36) peut s’écrire 
alors, en tenant compte de 2 et (36), sous la 
forme suivante : 
Div(QuW) 
> TO > 
. )/1+ € = D(Q, Q'— Q)d0 


AS 0: Jy/1+ . r(Q)+Q( (Q)]. (37) 


Or, comme on considère un courant des particules 


s 

en diffusion à vitesse Q, donc la vitesse des parti- 
cules considérée, dans le volume dV au moment f{ 
est supposée constante, on doit encore avoir 


Qu Grad W 


Pre RATE 
_ ASS Q, )y/i+ = [F(Q) + a(Q)] 


FO HART Die 
+] V Fe. pi, 054) 80,0 Q4@, 68) 


ou encore, d’après (16) et (25) 


7e ŒOW 

1 

\ ec? ot 
RERO 

<p/it <= | 
e 


x ae Le(Q, 0 = 0)40. 


(0. grad W AT HT es Où ) 


r(Q)+a(0)1æ+ [mr 00) 
(39) 


(5) F. JÜTINER, Ann. der Physik, 4° série, 1911, 84, p. 859. 
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L'équation intégro-différentielle (38, 39) doit être avons 
la généralisation relativiste de la formule de 
O. Halpern, R. Lueneburg et O. Clark (°). L’équation a 
donnée plus haut peut être écrite aussi dans VAE 
l’espace %, y, 2, Lx, V,, V, de la manière suivante Le \C 
+ - 
en introduisant les nouvelles fonctions w,(r, v, t), 
o (v, v'), y (0), © (0), , 


0% te ) 
+ (v, gradwæ 
dt 22e : É d 


: D Pa r ù Le 
PRE as Re ie > V4 ( 
Re Ù SAT, c? S ( 92 \° 

| . 
ee EN c? 


= 
ee PER p, t) [Y(c)+ow(c)]. (40) 


or On peut en effet exprimer l’élément de volume de TE 
7 _ _ = 
| vitesse dQ en fonction de l'élément dv (1) VA ( 
Et x 4 
dQ = LT (41) 
| (ÿ-3) 
f C° OI 
e =c(?) ed (5 (die = = = 
Complément de calcul. — a. Démonstration de PE No 5 RÉ 


“ la formule (36). — En effet, on a à = RU RE ne 
+ le système de coordonnées étant choisi lié au corpus 
As Ca. us by cule en mouvement dans la direction Ox; il s'ensuit 


Divisant par c? et en ajoutant l'unité aux deux 
membres, on trouve dOx= 


5 +. “(9 GITE 

3 Q? ne ES » : | 1 

ue — = +, RCRRRRE ee o : Ft 
de ee. ia Fe) - POSE Fer res le 5 des =: 

| A 1 (} ; p2 \2 € F. 
ie a. JE 


b. Démonstration de la formule (41). — Nous 4e : 


ie 

(°) O. HALPERN, R. LUENEBURG et O. CLARK, Phys. Rev. Chat | TS dre 
1938, 53, p. 173. | ( Ur 
ce 


(°) Pour la démonstration de cette relation comparer : 
J. BECQUEREL, loc. cit., p. 100-101; voir aussi F. JÜTINER 
Ann. der Physik, 4° série, 1911, 34, p.. 862. 


, 


Manuscrit reçu le 10 février 1940. . 
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